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1. Introduction 

The discovery of cosmic acceleration requires a source of new physics which may come as a modifica- 
tion to gravity, but is entirely consistent with a cosmological constant or quintessence, a dynamical 
dark energy [1, 2, 3].''^ Current observations suggest that dark energy comprises approximately 
70% of the energy content of the universe [7, 8] . 

^ Other suggested alternatives include the backreaction of inhomogeneities on the global expansion rate [4]; recent 
implementations of this have individually produced a small dark energy [5] and a large variance in the Hubble rate 
[6] but not provided a unique dark energy source. 
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The second major ingredient of modern cosmology is the dark matter, making up approxi- 
mately 25% of the universe's mass [7, 8]. Yet, remarkably, the existence of both of these dominant 
components of the universe is only deduced from their gravitational interactions. While a funda- 
mental microscopic picture is still lacking, a unified description of the dark sector might exist. 

One way to test such a proposition is to explore the phenomenological consequences of possible 
further non-gravitational interactions between a dynamical dark energy component, an ultra-light 
quintessence field, and dark matter [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]. As an immediate 
consequence of the interaction, the scalar field mediates a long-range force between the dark matter 
particles. As a result, for some scenarios strong constraints on the strength (3 of this dark force have 
been derived from dark matter dynamics [20, 21] and from the modified evolution of cosmological 
perturbations and subsequent impact on the CMB [22], constraining the coupling to 5% — 20% of 
gravitational strength. 

On the other hand, there are several observational motivations for dark interactions beyond 
gravity which challenge the ACDM concordance model. These include anomalously large bulk 
flows on scales of 100/i~-^Mpc and beyond [23], and higher supercluster densities [24] and sparser 
voids [18] than expected. 

Certain types of adiabatic coupled dark energy scenarios [25, 26, 27] have been shown to 
naturally allow for strong couplings /3 S> 1 [28, 29, 30, 31] between a light scalar and all forms of 
matter - even baryons - by means of the so-called thin-shell or chameleon mechanism. This applies 
if the field is kept in a minimum of its effective potential inside an overdense region and implies 
that deviations from the gravitational force law are shielded in the interior of bound objects. 

Another alternative for scenarios with couplings beyond gravity recently studied in the liter- 
ature is coupled neutrino quintessence [32, 33, 34, 35, 36]. The scalar /neutrino interaction results 
in growing neutrino masses which can trigger the onset of cosmic acceleration [36]. 

Since all viable dark energy models must yield extremely similar background behaviour, there 
has been intense study into alternative methods to discriminate between individual homogeneous 
and inhomogeneous models at a perturbative level. An uncoupled quintessence field with negligible 
self-interactions is not sensitive to the presence of inhomogeneities in the universe, and thus remains 
almost perfectly homogeneous. However, within overdensities the presence of a scalar/matter 
coupling causes the field to collapse along with the dark matter which itself feels an additional fifth 
force. The resulting modified growth of perturbations has been analysed in the linear [37, 38, 39, 40] 
and in the non-linear [41, 42, 43, 44, 45, 46, 47, 48, 49, 50] regimes. 

However, during structure formation, the dark energy provides a subdominant component to 
the energy budget of the universe, thereby limiting the size of the predicted effects. At late times, 
not only does the dark energy dominate the universe, but the matter distribution is significantly 
non-linear. It then seems natural to expect further sizeable effects at a non-perturbative level, 
emerging from non-standard scalar /matter interactions in the recent universe. In this paper, we 
therefore focus on possible corrections to the homogeneous equation of state of a coupled scalar 
field resulting from its non-trivial spatial profile acquired within clumped matter. 

As a result of its formation and collapse, a gravitationally bound structure in the universe 
is characterised by a non-trivial density profile with a shape which reflects its evolutionary state. 
By the current epoch, galaxies and galaxy clusters have had time to develop into fully virialised 
structures and, according to both observations and simulations, are well-described by a Navarro- 
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Frenk- White density profile. Even though these objects have a sizeable abundance, their total 
volume makes up a relatively small fraction of the Hubble volume. While it is certainly conceivable 
that local fluctuations in the coupled quintessence field could produce noticeable corrections to 
the homogeneous equation of state from such objects, we might expect the greatest corrections 
to result from the largest and most massive known structures in the universe, the superclusters. 
However, in contrast to fully relaxed objects on smaller scales, superclusters are still currently in 
the process of forming. This implies that the spherically averaged density profile of a superclustcr 
depends on its evolutionary state. In the future, it will tend for all superclusters towards an 
isothermal form proportional to at all radii [51]. However, at present the matter profiles of 
less evolved superclusters are typically shallower at small radii, but exhibit a steeper slope at large 
radii [52, 51]. 

Our first aim in this paper is to analytically solve the scalar field equations for different kinds of 
field potential and sourced by matter comprising of isolated, compact objects with density profiles 
of Navarro-Frenk- White and of isothermal form, corresponding to galaxies and clusters, and to 
superclusters respectively. The resulting spatial profiles of the field will allow us to determine the 
average equation of state of the scalar field within such bound structures. 

Our second aim is to carefully consider the impact of the evolutionary state of superclusters 
on the average equation of state of a coupled quintessence field. For this purpose, we will employ 
the results of numerical simulations to consider two case studies in different stages of formation. 
We will compare our numerical findings with those for an analytical fit to the simulated data, 
namely a Gaussian profile in the core stitched to an extended isothermal tail. 

Taken across the Hubble volume, these results will allow us to estimate the overall effect on 
the global equation of state resulting from the induced inhomogeneous nature of the coupled scalar 
field. 

After setting the stage for our analysis in §2, we define the average equation of state of an 
inhomogeneous scalar field. In §3.1 we provide a model-independent qualitative discussion of 
the scalar field dynamics in the presence of a coupling to clustered matter. In §4 and §5.1 we 
analytically determine the spatial profile of a scalar field coupled to bound matter structures 
with an NFW and an isothermal one profile, respectively, and calculate its average equation 
of state within such objects. In §5.2 we construct an analytical model for simulated density 
profiles of representative superclusters in different evolutionary states. In §5.3 we will discuss and 
compare our analytical and numerical results and estimate the implications on the equation of 
state averaged across the Hubble volume. §6 contains an application of our results to coupled 
quintessence scenarios with interactions beyond gravity and in §7 we conclude. 

2. Setting the Stage 

Our aim is to investigate the dependence of the average equation of state of a quintessence field (j) 
on its coupling to clustered matter. 

2.1 The Matter Distribution 

In the literature, various viable dark energy scenarios have been considered in which a light scalar 
field couples either to dark matter [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19], to neutrinos [32, 33, 34, 
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35, 36] or democratically to all forms of matter [25, 26, 27]. In general, all of these types of coupled 
matter have by the present time developed into bound compact objects and were shown to exhibit 
similar characteristic density profiles [46, 53, 50]. In order to deduce results valid for different kinds 
of coupled quintessence scenarios, we will therefore model the matter distribution as individual, 
isolated, spherically-symmetric clumps of a variety of sizes embedded in a background vacuum; we 
will consider bound objects of galactic size, as well as clusters of galaxies and super clusters. This 
is in contrast to the usual approach where overdensities are embedded in a constant background 
density. This implies that we will neglect non-linear interactions between the objects; since the 
typical distance between each of these bound objects is typically significantly larger than their 
size, this can be justified as a first approximation. 

We are likewise taking spherical symmetry to hold for objects across a wide range of scales; for 
virialiscd objects up to cluster scales we will assume a Navarro- Frcnk- White (NFW) profile [54], 
while for supcrclustcrs we investigate the characteristic isothermal profile along with modifications 
to it at small and large radii (see [52, 51] and references therein). Spherical symmetry is a good 
approximation on smaller scales where objects are fully virialised. When averaging across a large 
number of objects, we can assume it can also be justified on larger scales, and we would expect 
our calculations to provide reliable order-of-magnitude estimates. 

For the purposes of this paper it suffices to approximate the geometry of space as Minkowski, 
9^lv = 'Hfjii' — diag( — 1, 1, 1, 1), implying that both the Newtonian potential and the backrcaction 
due to the energy density in (p stay small everywhere. Even though we will be taking averages 
across the Hubble volume, this approximation is valid since we take the bound objects to be 
isolated from one-another; each resulting average is therefore taken only across a patch of space 
small compared to the Hubble volume. 

In addition, we assume the field dynamics to be static, since the characteristic time scale for 
a dark energy candidate is of the order of the Hubble time. Furthermore, in order to probe the 
dependence of the local profile of on the fundamental potential of the scalar field, we will consider 
two fiducial examples, an inverse-power law and an exponential potential. 

2.2 The Average Equation of State 

In this section, we isolate the relevant quantities that define the average equation of state of a 
scalar field with a self- interaction potential V{(j)). Since an interaction with clustered matter 
will induce local variations in the field value, we do not a priori assume that spatial gradients 
vanish under averaging as they do in the usual homogeneous approximation. 

When considering an "average equation of state" , there are three definitions we can readily 
imagine. The equation of state of a scalar field is defined from its pressure and energy density, 
which themselves depend directly on the scalar field itself. We can therefore construct an average 
equation of state wi by forming the local equation of state and averaging this across a spatial 
volume; or we can average the scalar field ^ itself and define the average equation of state W2 from 
the general forms of the pressure and energy density; or we can average the pressure and energy 
density independently and define an average equation of state from these. The three definitions 
can be written 




(2.1) 
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and they coincide for a homogeneous field. Since the relevant thermodynamic quantities are the 
pressure and the energy density rather than the scalar field or the local equation of state, we choose 
the last approach. 

Varying the scalar field Lagrangian with respect to the metric and averaging over the Hubble 
volume Vh, the average pressure and energy density of the field can be defined as 



P4> 



d^x 



2^ 



-A2 . 



+ 2 i^^r + ^i^) 



l-TT 1 



In the static limit ^ 2± 0, the average equation of state w therefore takes the form 



w 



where 



_ Pi 
Qct> 



-i-H_g 



(2.2) 



(2.3) 



(2.4) 



Accordingly, we can already isolate two exact limiting cases: for Q S> 1 the scalar field dynamics 
will be gradient dominated, leading to w — )• —1/3; conversely, for Q ^ 1 the field gradient is 
subdominant compared to the potential, and w ^ —1. This implies, therefore, that in the static 
limit -1 < ItJ < -1/3. 

Gradients of the scalar field originating from an interaction with clustered matter can thus, in 
principle, lead to sizeable corrections to w as compared to its uncoupled, homogeneous analogue 

Wh, 



_ P0 
^h = — 7=r 

Q<l> W 



-V{<f>) 



(2.5) 



where again the last equality applies in the static limit. 

To transfer the calculations for the equation of state of the coupled field within a bound object 
to an average across the Hubble volume it is necessary to make some assumptions concerning the 
distribution of matter sources. We will separate the bound objects into distinct classes, each of 
which produces a particular scalar field profile and therefore the same equation of state. If i 
denotes a particular class of objects with radius r-Qbj,?, then Eq. (2.3) can be approximated by 



W 2±= 



(2.6) 
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Qi 
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Here Ni is the number of objects in the class i, we have assumed Q — in the intergalactic medium, 
and angle brackets refer to quantities averaged over the volume of an isolated object Vobj.i- Note 
that TZi is a measure for the volume dilution of the contribution Qi to the sum. We can then 
immediately deduce that notable differences between Eq. (2.7) and Eq. (2.5) could well emerge 
from a class of objects for which the gradient contribution Qi is large enough to compensate for 
the volume dilution implied in TZi. 

3. Dynamics of the Scalar Field 
3.1 The Equation of Motion 

As is clear from the previous section, to estimate w we must determine the spatial profile of the 
scalar field (f> resulting from its interaction with clustered matter. While we assume standard dark 

matter density distributions for the coupled matter species, wc will consider different possible 
forms of coupling functions leading to distinct dynamics of the field. By this means our results can 
be generalised to hold for the various different kinds of coupled scalar field scenarios mentioned in 
§2. 

In general, in the presence of an interaction of strength /3a between ^ and a > 1 non-relativistic 
matter species, the dynamics of (f) are governed by an effective potential of the form 

yes{^,r) = V{<j>,r) + ^Pa{cP,r), (3.1) 

a 

where the first term arises from self-interactions of the quintessence field, generating a monotonically- 
decreasing potential t^(0) of a runaway form. The second contribution to V^ff originates from the 
<p dependence of the masses of the coupled species, ma{4>), introduced by the interaction. Since 
the corresponding energy density pa is a function of the particle mass, as a direct consequence 
Pai'iTT'ai^')) Contributes to the effective potential for 0. 

For simplicity, we assume a universal coupling of strength P = Pa and an effective energy 
density p leading to an effective potential of the form 

yeff(0,r) = F(</.,r)+p(0,r), (3.2) 

which encodes the effects of one or more coupled species on the scalar field potential. 

Note that depending on the sign of the coupling function dlog(m((;l!)))/d0, we must in general 
discriminate between two cases which lead to distinct dynamics. Namely, a positive sign allows 
the effective potential to exhibit a minimum, while a negative sign implies a (monotonic) decrease 
of VeS with (p. 

We take for simplicity the mass of the coupled matter to depend only very weakly on changes 
in the scalar field value such that it can be taken to be approximately constant even within 
environments with strongly varying matter density - that is, p{(t>,r) 2± p{r). Under the further 
assumption of spherical symmetry for the matter profile, the dynamics of the field in the static 
limit are described by the equation of motion 

0"(r) + -0'(r) = V^{4>, r) + /3p(r) = V,s,^{ct>, r) , (3.3) 
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where here and in the foUowing primes denote derivatives with respect to r, a subscript , </> is a 
derivative with respect to (j), and the last equahty results from Eq. (3.2) and a coupling function 
of exponential form. 

It will transpire that our results are not sensitive to the form of the potential V{(t)). We 
demonstrate this by employing two fiducial examples. The first is of exponential form, 

= M^ie"^'^, (3.4) 

where Mpi denotes the reduced Planck mass and a > 10 to fit observation [55]. 
Our second fiducial potential is an inverse power law potential of the form 

^(^) = ' (3-5) 

where n is a positive constant and M an intrinsic mass scale. It fixes the behaviour of ^ in 
the intergalactic medium and has to be compared to the scale of dark energy in a homogeneous 
universe. 

However, unless specified the potential should be taken to be general. 
3.2 Qualitative Description of the Solutions 

Let us start our analysis with a model-independent description of different possible kinds of ap- 
proximate analytical solutions to Eq. (3.3). 

We first describe the initial conditions that we assume to hold for any effective potential. 
Since the field equation in Eq. (3.3) is a second order differential equation, we have to specify two 
boundary conditions for (f). Firstly, to guarantee a non-singular behaviour of (j) at the origin, we 
require as usual 

= at r = . (3.6) 

To set the second boundary condition, we note that in our setting the bound matter clumps are 
embedded in a vacuum, and so the matter density outside of these asymptotes to zero and not to a 
non-zero background value as it would for an idealised homogeneous universe. We must accordingly 
take the field to be essentially free in the intergalactic medium, that \s (f) = (f){ree- Therefore, in 
'empty space' its dynamics are solely governed by its self-interaction potential V{(f){ree)- Our second 
boundary condition is then 

(j) — >■ 0free as r ^ oo . (3-7) 

To a good approximation, the solutions for (j) turn out to be independent of the scalar field 
potential V{(j)). However, there is one situation in which this is not the case, where (p — 0min; 
with ^min being the field value at the minimum. This can occur when the following conditions are 
satisfied: 

VeS,(p = 0, such that Vcs exhibits a minimum (3.8) 
i(0) <C <?!>min(0), where (f)i is the value of cj) at the origin (3.9) 



J 



P 



<^ m^{r), which we refer to as the adiabaticity condition. (3.10) 
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If these conditions hold, then the equations of motions ahow for an adiabatic solution to Eq. (3.3), 
valid in the range < r < rkick for rkick < ^obj, where the field instantaneously adjusts to spatial 
changes of the minimum Ves^ff,. Here, rkick denotes the radial coordinate at which the fields gets 
kicked out of the minimum. The adiabaticity condition corresponds to the requirement that the 
matter density changes on scales much larger than the Compton wavelength of the scalar 
field in that region. 

Since in this case the source term is initially V^g^^ — 0, and the mass of the field inside the 
object is large, m(^(0) » rrif^cc, perturbations in the field will be exponentially suppressed. As a 
result, if r < rkick and Eq. (3.10) is valid, the matter source term p(r) can be treated as a constant 
in the scalar field equation of motion [26] . 

As soon as r > r^ck, condition (3.10) is violated and the field gets displaced from the minimum, 
and the contribution from to the field equation, Eq. (3.3), becomes negligible with respect to 
the matter density p, 



Ml 



PI 



(3.11) 



Note that if the condition (3.9) is violated, this implies that the field is already sufficiently 
displaced from the minimum at the origin such that rkick = and its evolution is non-adiabatic 
from the start. 

As it will turn out, for the density profiles and potentials discussed in this paper, Eq. (3.10) is 
already violated for very small r. Accordingly, the field is immediately kicked out of the minimum of 
Vcg even in scenarios that in principle might obey all three conditions. As a result, the evolution 
is effectively non-adiabatic across virtually the entire object. Our result agrees with previous 
numerical studies, which found that the so-called chameleon effect [25, 26, 27] is not significant for 
compact structures at low redshifts [56, 57, 58]. 

In general, therefore, in the scenarios we will consider we are lead to the second type of 
approximate solution to Eq. (3.3). This solution is valid as long as Eq. (3.11) is fulfilled and 
implies that Eq. (3.3) reduces to 

(3.12) 

r Mpi 

Following the approach described in [28], we have verified that this approximation stays valid for 
< r < robj in the scenarios under investigation in this paper. 

Note for later reference that Eqs. (3.6) and (3.12) yield the following result for the gradient 
squared of <p, 

(V</>(r))2 = [j^^'^'' drfV(r)j , (3.13) 

independent of the form of the scalar potential V((p). 

Averaged over the volume Vobj = ^'^''^ohj/^ object, this gives 

2 / \ ^ 

((V0)2) = AttV;^] (^^) I drr-2 J dffV(r) • (3.14) 
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Finally, outside of the object the dynamics of the field are only governed by its self-interaction 
potential V with <C 1. Accordingly, in this regime the field equation reduces to 

(^" + -(/)'~0, (3.15) 
r 

which gives the solution 

Ce(~"^free(r-r-obj)) 

0(r) = (/(free • (3.16) 

r 

Here we have employed that the Compton wavelength of the field in the intergalactic medium, 
m^gg, is much larger than the radius, Tobj, of the object. Furthermore, the constant C can be 
determined by requiring continuity of (p and d^/dr at r = robj- 

4. NFW matter density profile 

We now turn to objects of galactic and cluster scales, for which we employ an NFW matter density 
profile [54] , which is a phenomenological fit to the density distribution of dark matter halos in CDM 
cosmological simulations, 

with 

_A c^pc _ rvir 

3 1og(l + c)-c/(l + c)' r, ' ^ ^ ^ 

and we have defined the dimensionless variable R = r/rg, where < i? < c, with r being the radial 
coordinate and the characteristic scale of the halo. Together with its virial radius ryir, Tg defines 
the concentration parameter c. Furthermore, the characteristic density ps is proportional to the 
the critical density of the universe, Pc, and to the virial overdensity. A, which can be computed 
as described in [59]. While the two parameters ps and are in principle independent, they where 
shown in simulations [60] to be correlated through the concentration parameter c, the critical 
density pc and the virial mass Myir, 



Mvir = 4Trpsr, 



log(l + c) 



1 + c 



(4.3) 



Note that typical virial masses of galaxies and galaxy clusters are Myir — 10^" — IO^^Mq and 
Mvir ^ 10^-^ - W^^Mq, respectively. 

In the following, we will employ the relations [61] 



As discussed in the last section, the behaviour of the coupled scalar field inside a matter lump 
crucially depends on the magnitude of its interaction range m^^ compared to the rate of change 
of the matter density, see Eq. (3.10). 

In the following subsection we derive the range of validity of the adiabatic approximation to 
the field equation for the two fiducial potentials in Eq. (3.4) and Eq. (3.5). In addition, we assume 
the matter coupling function to be an increasing function of the field value such that the effective 
potential exhibits a minimum. 
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4.1 Adiabaticity Condition for an NFW Matter Profile 

For the NFW distribution the adiabaticity condition (3.10) becomes 



p'{R) 



p{R) 



1 l + 3i? 1 _i 

~ —R ^ , (4.5) 



Ts R{1 + R) rs 



where the last equaUty assumes i? <C 1. As a result of the cusped structure of the matter density, 
Eq. (4.5) diverges at the origin. This already sTiggests that the field will be immediately displaced 
from the minimum, as we now explicitly verify. 

In the case of the exponential potential, we find small excitations about the minimum of the 
potential to be given by 



1 



m,iR)=^-^il + Rr^R-L (4.6) 



The resulting adiabaticity condition therefore takes the form 

12Mpi \5l ^ ^ 

1+11- -^—^ 1 , (4.7) 



-Tt ^ — — + 



3 18 V Mpi 

with 



aA^NFW 



^*-w.^ = O(:0-)M„(^)°", (4.8) 

where the last equality results from Eqs. (4.3)-(4.4) and holds for c > 1. As it will turn out later 
on, A(/)NFW is a measure for the absolute change of the scalar field value within the object, which 
for typical model parameters is much smaller than the characteristic mass scale in the potential, 
A^NFW ^ Mpi/a. Eq. (4.7) therefore cannot be solved for any R> and the field evolution is 
non-adiabatic from the start. 

In the case of the inverse-power law potential, the field mass expanded around the origin takes 
the form 

1(11+2) 

m,{R) c V^n-K^)M-i(^) (^^) " • (4.9) 

Again, for any positive n, Eq. (4.9) diverges less fast than Eq. (4.5) at small radii. 
In this case, the adiabatic condition in Eq. (3.10) close to the origin becomes 

, t n + l \ 1 n + 4 —2 /n-K^^ 

< (ra + l)'^—)n~nM~rs " (A^nfw)"*^"^ < 1 , (4.10) 

where the last equality holds even for the largest and most massive clusters of galaxies and realistic 
model parameters. 

We therefore conclude that Eq. (3.11) can be assumed to hold and the field evolution is non- 
adiabatic within objects up to cluster scales at the current epoch. This result is in agreement 
with corresponding numerical studies that have demonstrated the absence of the thin-shell effect 
for compact structures at late times. In a series a papers [56, 57, 58] the authors employed full 
numerical simulations to analyse the development as well as the disappearance of thin shells in 
chameleon-like scalar field scenarios as a function of redshift during nonlinear structure formation. 

Note also that the resulting solutions for the field thereby become independent of the functional 
form of the scalar field potential. 
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4.2 Non-Adiabatic Solution for the NFW Matter Profile 
For an NFW profile, we find that the solutions to the field equation in Eq. (3.12) are 

'log(l + R) 1 



4>m. = <j>h:ee — A^nFW 



R 



1 + c 



^out 



^free - A^NFW log(l + c) 



(4.11) 
(4.12) 



l + cj R 

tree'vir ^ 1- Furthermore, as described in §3.2, we have deter- 



where the result again assumes mfree?^v 
mined the integration constants by matching the inner solution to the outer solution ^out) 
R = c and required continuity of the first derivative of the field at i? = c. 

As we mentioned in the last section, the absolute change in the scalar field value, (f){c) — (p{0) ~ 
A(/)NFW defined in Eq. (4.8) within a compact object with virial radius Tvir, is a measure for the 
total mass Mvir defined in Eq. (4.3), which is interacting with with coupling strength /3. 

Using Eq. (4.11), the volume average of the exponential potential becomes 



{V) ~ Ffr, 



1 + 



NFW 



Mpi 



\ + {l-c) (2-31og2) 



(4.13) 



for c > 1. 

In the case of the inverse power law potential we find 



1 + n 



A^NFW f 1 



+ (l-c) (2-31og2) 



(4.14) 



for c > 1. Note that for both potentials, the matter coupling results in a positive correction term 
with respect to the free field potential T4rce- However, importantly, for typical model parameters 
applying to galaxies and galaxy clusters, it turns out to be negligible, since according to Eq. (4.8) 
A0NFW ^ 0free a-^d A^nfw ^ Mpi/a, respectively. 

For the volume average of the gradient squared we get 



{iVcf>f) = 3c 



-3 f A^NFwV 



dR 



log(l + R) 
R 



1 



1 + R 



3/(c 



/ A^NFW Y 

V crs J 



with 



/(c) = 



1 + c 



log(l + c) 



(4.15) 



(4.16) 



where /(c) ~ O(10~^) for typical values of the concentration parameter c. Note that even though 
the matter density formally diverges at the origin, the gradient squared stays well-defined dX R = 0. 
This is due to fact that the volume integral over the matter density is finite as can be seen from 
Eq. (3.14). 

Accordingly, using Eq. (4.4), the volume average of the gradient squared normalised to the 
free field potential becomes 

\2\ / )ij \ 0.68 



Q 



((V'/')0 

^free 



5.12 X 10-^/3^ 



(4.17) 
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for typical galaxies, and 

for higher mass clusters, such that Q ^ 1 for couplings of (sub-)gravitational strength ,0 < 1 (see, 
however, §6, where wc analyse the results for coupled quintessence scenarios with interactions of 
super-gravitational strength). 

Recall that, according to Eq. (2.7), contributions to W depend on the product of Q and TZ. 
However, for all classes of objects on galactic and cluster scales, the volume suppression amounts 
to at least TZ < O{10^^^)N <^ 1 for typical abundances of these object classes, w therefore receives 
negligible corrections compared to the uncoupled case Wh — —1- 

As a concluding remark, let us note that the small result for Q is a direct consequence of the 
specific functional form of the NFW profile. Namely, according to Eq. (4.15), the integrand of 
the volume average is of 0{R^) for R 1 and of 0{1/R^) for R^ 1. This implies that in both 
regimes the corresponding contributions to the volume average are strongly suppressed and thus 
the average gradient squared only receives sizeable contributions from radii i? ~ 1 <C c. 

However, as we will see in the next section, for the density distribution characteristic for 
objects on supercluster scales, the average field gradient receives almost equal contributions at all 
radii and thus turns out the be much larger. 



5. Superclusters 

In this section we investigate the impact of superclusters, the largest known bound structures in 
the universe, on w in coupled quintessence. Within a ACDM cosmology, superclusters are expected 
to have masses ranging from Mgc ~ lO^^/i~^M0 up to greater than 6 x 10^^/i~^MQ, with sizes 
~ 10/i~^-160/i~^Mpc [52, 62]. Superclusters therefore fill approximately 3% — 5% of the current 
Hubble volume [52]. 

Since structure forms hierarchically, the superclusters we observe in the present universe are 
still in the process of formation, while galaxies and clusters of galaxies have already developed into 
collapsed and virialised structures. Accordingly, depending on their evolutionary state, superclus- 
ters tend to have a much shallower density profile than fully relaxed objects. 

According to simulations, the spherically averaged density profile of superclusters is well-fitted 
by an isothermal profile p{r) oc over a broad range in radius. While the profile is shallower 
than this at low radii, it steepens at larger radii [52, 51]. Since we only have statistical information 
on the exact density distribution, we will proceed in the following manner. 

To establish some intuitive understanding of the behaviour of the field inside a typical super- 
cluster, we will begin by analytically determining the generic profile of a scalar field interacting 
with isothermally distributed matter. This allows us to get a feeling for the magnitude of the 
effects on the global equation of state as a function of the coupling, the mass and the radius of a 
single, standard supercluster. 

Subsequently, we will turn to two simulated supercluster case studies. These allow us to test 
both analytically and numerically the dependence of our results on deviations from the isothermal 
profile for small and large supercluster radii. 



- 12 - 



Finally, we will end our discussion of superclusters by estimating the overall effect on the 
global equation of state. 



5.1 Isothermal Matter Distribution 

Formally, the isothermal profile p{R) oc exhibits a singularity at the origin, where R = r/rsc 
denotes the radial coordinate normalised to the supercluster radius rsc- However, bound objects 
which have not yet relaxed will exhibit a flattening of the matter distribution at some R = A with 
< A <C 1. Accordingly, for A < R < 1 the isothermal distribution expressed in terms of the 
finite matter density p = p{A)A'^ reads 



With the help of Eq. (3.11), the field equation Eq. (3.3) then takes the form 

d'^d) . 2 d</> /3 



di?2 RdR 



P_ 

Mpi i?2 



Assuming both the field value and its first derivative to be continuous at R 
A< R<1 that 



4>m = ^free - A0is 
^free - A0is 



1 - log(i?) - 



A 
R 



R=A/ J 



1-A I 



where 



Mpi 



R=A 



(5.1) 
(5.2) 

A, we get for 

(5.3) 
(5.4) 

(5.5) 



As described in §3.1, we have determined the integration constants by matching the inner solution 
(l)in to the outer solution (j)out, and by requiring continuity of the respective derivatives at i? = 1. 

Note that we can disentangle the different contributions to the average gradient squared by 
investigating the limiting case A0is S> A d(?!>/di?|^^^. Doing so assumes the profile of (j) in the 
isothermal region to be insensitive to its slope in the core region; in general, this slope is associated 
with large uncertainties introduced by the matter profile in the core. Comparing this with the 
analysis in §5.2, where we employ a Gaussian distribution at small radii, will allow us to gain some 
understanding of the influence of the core profile. 

According to Eq. (5.3), for A(^is S> A d(j)/dR\j^^y^ we find the gradient squared of an object 
with isothermal density profile averaged over the volume of the object to be 



{iVct>f) = 3 



IS 



rsc 



1 



/ 



dR(^l 



A 
R 



J AMY 
V ^sc / 



1 + 



21og(^) 
A 



^ 3 (^) ' [1 + 2 log{A)A + 0(^2)1 , 

V ^sc / 



(5.6) 



where the last equality assumes A <^ 1. Note that in this limit the integrand of the volume average 
is of order unity for ^ <C i? < 1 implying that the average gradient squared receives almost equally 
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large contribution at all radii in this range. Contrasting this result with its analogue in the NFW 
case in Eq. (4.15), we observe that the isothermal matter distribution results in an average gradient 
squared which is a factor /^^(c) ~ 0(10^) larger. We therefore arrive at the crucial result that the 
impact of the interaction on the coupled scalar field equation of state strongly depends on the radial 
matter distribution. Since the spherically averaged supercluster profile is typically isothermal over 
a broad range in radius, this result already suggests that the dominant contribution to the average 
gradient squared of the field will originate from this region. 

5.2 Composite Supercluster Profiles 

In the quantitative investigation that follows, we consider two representative superclusters, SC 8 
and SC 98, with masses of Msc ~ 5.4 x lO^^h'^MQ and Msc ~ 3.6 x lO^^/j-^Mo, respectively, 
taken from a simulation performed on the Beowulf Cluster at the University of Groningen [51]. 
While at the current epoch, SC 8 is already a centrally concentrated object, SC 98 is still more 
extended. 

However, both supercluster density profiles exhibit as a characteristic basic shape a high- 
density core with a central density of ~ (10*^ — I0^)pc, embedded within an isothermal power law 
region with slope —2. 

Furthermore, for both of these examples we find that the central region is well-fitted by a 
Gaussian density distribution. Since the two examples represent in some respect two extremes of 
supercluster evolution, it is a reasonable assumption that this is a generic feature. In this section, 
therefore, we model a generic supercluster density profile by stitching together at low radii a 
Gaussian and an isothermal matter density profile. This then allows us to determine analytical 
solutions for the scalar field profile within typical high-mass superclusters which will be compared 
to our corresponding numerical results for the spherically averaged density profiles of SC 8 and 
SC 98. 



5.2.1 Gaussian Matter Distribution 

In the following we determine the field profile (pQ arising from a small Gaussian matter region 
stitched to an extended isothermal tail. The Gaussian and isothermal regions are connected at a 
small radius yl <C 1. 

We parameterise the Gaussian density profile by 

PG(i?)=Pmaxe-^(^^''^% (5.7) 

with 



A 



(5.8) 



We have defined the dimensionless constant Cg = robj/c, a being the half- width of the Gaussian 
distribution with amplitude A. 

In the following subsection we derive the range of validity of the adiabatic approximation to 
the field equation for the two fiducial potentials in Eqs. (3.4-3.5) and will argue that the results 
hold for any realistic analytical scalar field potential. 
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5.2.2 Adiabaticity Condition for a Gaussian Core Profile 

If the effective potential exhibits a minimum and l^(^) is of exponential form, small perturbations 
about the minimum are described by 



Accordingly, the minimum condition in Eq. (3.10) reads 

Cg Mpi ' 

where the last equalities hold for typical parameter values for which aA(j)G/Mpi <C 1 and Cg S> 1. 
Here W{x) is the Lambert function, which satisfies W{x)e^^^^ = x. 

If the potential exhibits a minimum and V{(p) is of inverse power law form, we get 

m^{R) = Ic(^)e-'^('^^y^^^\ (5.11) 

with 

JC = x/^hH: (nM"+^) ±^ A(Pl}"+'^ . (5.12) 

\fobj J 

With the help of Eq. (5.7) the minimum condition in Eq. (3.10) becomes 

(5.13) 

Furthermore, the last equalities hold for /C ^ 1 and Cg 1 which will turn out to be a good 
approximation for typical parameter values (see §5.3). 

As we will demonstrate explicitly in §5.3, the minimum solution in practice has a very limited 
range of validity, since the evolution turns out the be non-adiabatic very close to the origin. 

5.2.3 Non-Adiabatic Solution for the Gaussian Core Profile 

In the non-adiabatic case, in the core region, < i? < ^, the field equation takes the form 

d^^2#^^^ 

with 



-^e^l-^'.*' , (5.15) 

where the last equality results from matching Eq. (5.1) and Eq. (5.7) R = A. 

Requiring d(?!)/di?|^^Q = 0, with the help of Eq. (5.3) and Eq. (5.5) for < i? < ^ we arrive 

at 



-{1 + (Ac, flo,i A)) e-i'-'> f 



(5.16) 



^'9 

AcPc = ^(-] =^^-^A0is, (5.17) 

Mpi \CgJ {ACgY 

(5.18) 
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where the relation in Eq. (5.18) guarantees eontinuity of (j) and its derivative at i? = ^. 

We now have an analytical solution for which allows us to determine the slope of the 
isothermal profile at R = A. 

5.2.4 Isothermal Distributions with a Gaussian Core Region 

In this section, we determine the dependence of ((V(/))^) in the isothermal region on the slope of 
the field profile in the core region. For this purpose, we require the scalar field profile and its 
derivative to be continuous at the matching point A. Recall that this condition fixes the relation 
in Eq. (5.18). 

Taking the derivative of the field profile in Eq. (5.16) at R = A and inserting it for d^/di?!^^^ 
in Eq. (5.3), we arrive at 



(t^in = 0free - ^^IS 



log(i?) - -X 



(5.19) 



0free - A(/)is [1 - Ax] 



-mfrcorsciR-^) 

R 



where 



X= 1 + 



erf e5(^-«)' 



The scalar field profile in Eq. (5.19) leads to a gradient squared averaged over the volume of the 
object of 

2 r 



((V<^)') 



IS 



rsc 



i + ,?M!)_i)^V + (x=-i)^ 



(A) - - ]a + o{a'^ 



(5.20) 



where in the last equality we have employed Eq. (5.18) and assumed A <C 1. Noting that for 
X — 7- 1 we recover the result in Eq. (5.3) with A d(/)/di2|^^^ <^ A4>is , we are now in a position to 
estimate the impact of the field derivative at the matching point on the average gradient squared. 
For ^ ^ 1 and log(yl) <C 0, the field gradient at R = A leads to an increase of the result by a 
factor ~ (1 — 2/3log{A)A) > which amounts to 0(15%) for typical parameter values. 

Accordingly, we find that a small central region with (approximately) Gaussian matter dis- 
tribution has a sizeable effect on the scalar field dynamics. This has to be contrasted with the 
negligible impact of the corresponding NEW matter distribution at small radii which was found 
to result from its cusped structure. 

For the volume average of the potential in the case of an exponential potential we get 



{V) ~ Ffr, 



1 + 



A 



aA(f> 



NFW 



Mpi 



In the case of an inverse power law we arrive at 



{V) ~ Ffr, 



4 \ A(/)NFW 
A \ 11 — 7 

3 / (pfree . 



(5.21) 



(5.22) 
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In both cases, the quintessence potential could be expanded around the free-field potential for 
aA0/Mpi <C 1 and A^/^free ^ 1, respectively. As in the case of the NFW profile in §4, we find a 
negligible correction term. 

Note that our result can be generalised to any analytic quintessence potential which allows for 
such an expansion. So long as A0 is much smaller than the characteristic mass scale of the scalar 
potential, the correction term will remain negligible. 

In summary, for all compact objects from galactic up to supercluster scales, the absolute 
change in the scalar field value induced by the scalar /matter coupling has a negligible effect 
on the average potential of the scalar field. Taken across the Hubble volume, this implies that 
possible corrections to the homogeneous average equation of state can directly be traced back to 
the square of non-trivial field gradients. Note also that as a result of the small present value of 
V — ^frcc — 10~^^'^Mp[, the field gradients only have to be of comparable magnitude in order to 
induce sizeable effects on the average equation of state. 

5.3 Supercluster Case Studies: Analytical and Numerical Results 

In this section we discuss our analytical and numerical results for SC 8 and SC 98, which we take 
as representatives for further and less evolved superclusters respectively. The difference in their 
evolutionary state is refiected in the shape of their spherically averaged density profiles, plotted in 
Fig. 1 together with a Gaussian (for i? <C 1) and an isothermal profile (for R < A) fitted to the 
data. The distribution of SC 8 drops faster than at i? ~ 0.5; in the case of SC 98 the slope 
is already changing at a much smaller radius, R ~ 0.1. Note also that SC 8 has a larger mass 
than SC 98 and therefore its impact on the the scalar field value is larger. This is demonstrated 
in Fig. 2, which shows the scalar field profile for the two superclusters according to Eq. (5.16) and 
Eqs. (5.19-5.20). 

We are now in a position to verify the limited range of validity of the adiabatic solution 
according to Eqs. (5. 10) and (5.13). For our test superclusters we arrive at 



for SC 98. Therefore, for both cases, the adiabaticity condition is already violated at i? ^ 1 such 
that in practise the field evolution is non-adiabatic. 

In the following, we will compare and discuss our analytical and numerical results for Q as 
defined in Eq. (2.7) and analyse the impact of deviations from the isothermal profile at large radii. 
Finally, based on our results, we will estimate w. 

Taking Fig. 1 as a rough guideline for the relevant scales, according to Eq. (5.19) we find for 
the isothermal matter distribution that 



A^G = 1.8 X IQ-^Mpi 



(5.23) 



for SC 8, and 



A(f>G = 5.8 X lO^^Mpi 



(5.24) 




(5.25) 
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Figure 1: Spherically averaged matter density profiles of the further and less far evolved superclusters 
8 and 98, respectively. Short dashed lines denote the Gaussian core region, long dashed lines the R^^ 
isothermal tail, and solid black lines the numerical results. 

To compare this result with our numerical findings, we insert the fitted radial density distri- 
butions shown in Fig. 1 into the general solution for the average gradient squared in Eq. (3.14)^. 
Fig. 3 shows our numerical and analytical results for the normalised average gradient squared, 



and its integrand of as a function of R. As demonstrated in the lower panels of Fig. 3, the final value 
of the normahsed average gradient squared for SC 8 is Qs = 1-63, almost an order of magnitude 
larger than that for SC 98, where Qgg = 0.21. Furthermore, for the former the analytical and 
numerical curves are in very good agreement for all radii, while in the latter case the analytical 
result overestimates Q{R) even at small radii. As can be understood from the plots in the upper 
panel of Fig. 3 of the integrand of Q in Eq. (5.26), the reason is that the isothermal profile overall 
provides a better fit to the profile of the further evolved SC 8. 

According to the analytical result discussed at the end of §5.1, Q receives almost equal con- 
tributions at all radii in the isothermal region, which is reflected by the small slope of the dashed 
curves. As could be expected from Fig. 1, the numerical curves agree with the analytic up to the ra- 
dius at which the numerical profiles drop faster than R~^. However, for larger radii, the numerical 

■^Note that we have extended the data at small radii with the analytical data from the Gaussian fit in order to 
employ the initial condition in Eq. (3.6) for at R = 



R 




(5.26) 







- 18 - 



m FT 

*free 1 " 

0.99 - 

0.98 - 

0.97 - 

0.96 




i 1 1 1111 




y / 


■ 


y / 




/ / 


_ 


/ / 




y^ 1 


■ 


■ — '"^ 1 


■ 


1 


- 


1 


■ 


/ 




/ 




/ 




/ 




/ 




y 




1 1 1 .... 


..,.1 . . 



01 



0.1 1 
R 



10 



Figure 2: The scalar field profile (/)(i?)/(/)fioc for the analytie fits for SC 8 (lower curve) and SC 98 (upper 
curve) and 13 — 1. The solutions in the Gaussian region are plotted with short dashed lines, those in the 
isothermal region with long dashed lines, and the external solutions with solid black lines. 

contributions to the integrand Q{R) show a decaying behaviour. Therefore, to a good approxima- 
tion the corresponding maximal value of Q{R) fixes the total Q, since further contributions are 
negligible. 



Based on our results in this section, we are finally in a position to estimate the overall effect on 
the expansion rate for coupled quintessence scenarios. To achieve this, we note that a particular 
class of high-mass superclusters presumably yields the dominant impact on the expansion rate; 
superclusters tend to become ever less concentrated the larger (and less evolved) they are, as 
reflected by their mean density, 

(5.27) 



3 Msc 
47r rl^ 



which implies that ((V(/))^) oc (/3/Mpi) M^^/r^^. Here, the factor of proportionality depends on the 
radial distribution of mass within the object. Adding to this, the abundance of bound structure 
drops with increasing radius of the object. We therefore conclude that the main source of scalar 
field gradients will stem from medium-sized superclusters which are both large and numerous 
enough that their volume still makes up a reasonable fraction of the Hubble volume (see for 
comparison the discussion in §2.2 below Eq. (2.7)). 



As a result, we assume the total volume taken up by all superclusters (~ (3 — 5)%Vh) to be 
mainly provided by one representative class of superclusters with a medium size of ~ 30/i-^Mpc 
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Figure 3: R^{V(f)f (top) and Q{R) (bottom) for SC 8 (left) and SC 98 (right) and /3 = 1. 

and normalised gradient squared Q. According to Eq. (2.7), this corresponds to TZi = TZ ~ 
(3 — 5) X 10^^. In Fig. 4 we plot the resulting w for two different values of 7^ as a function of Q 
together with the corresponding values for W which we get for superclusters similar to SC 8 and 
SC 98. 

From these plots we can see that, if superclusters similar to SC 98 dominate, only small devia- 
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Figure 4: w as a function of Q for a volume of 3% (left) and 5% (right) and /3 — 1. 



tions from uJ ~ — 1 are to be expected, even should the superclusters fill 5% of the Hubble volume. 
However, for further evolved superclusters of type SC 8, the deviations from the homogeneous 
equation of state are dramatic. For a lower abundance, where the superclusters fill 3% of the 
universe, the equation of state is pushed to IZ; ~ —0.98, slowing the cosmic acceleration to some 
degree. For the higher abundance of 5%, the equation of state is driven to ItJ ~ —0.97. Let us 
remark that the authors of [52] have performed a complete analysis of a dark matter simulation 
taking into account the shape, multiplicity, size, and radial structure of superclusters. They found 
that across the entire sample the spherically averaged density profile is roughly isothermal in the 
range 0.3 < -R < 0.7. This investigation suggests that on average, superclusters tend to be much 
further evolved than SC 98 and, in fact, exhibit a similar evolutionary state to SC 8. In light of 
these results, we therefore conclude that an overall correction of > 3% to the homogenous equation 
of state is conceivable in scenarios with scalar /matter couplings of gravitational strength. In the 
next section, we comment on the class of coupled scalar field scenarios which allow for (3 of 0(1). 
Furthermore, in §6 we provide an outlook of the implications of our results for scenarios with cou- 
plings of super-gravitational strength such as coupled neutrino quintessence [32, 33, 34, 35, 36]. In 
summary, as anticipated, in coupled quintessence scenarios the main source for notable corrections 
to w stems from the largest and most massive bound objects known. Since superclusters are typ- 
ically less than two orders of magnitude smaller than the Hubble radius, the volume suppression 
is many orders of magnitude lower than that for galaxies and galaxy clusters. 

In addition, the spatial changes in the coupled scalar field value were found to get more 
pronounced for higher masses of bound objects and are also sensitive to its distribution. Namely, 
for superclusters, the dominant isothermal region exhibits a shallower slope than the NFW profile 
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discussed in §4 and, importantly, lacks a very highly-concentrated central core (the cusp). Thus, 
according to their evolutionary state, the cumulative mass of an supercluster and its impact on the 
scalar field value change gradually over a much larger range in radius. This has to be contrasted 
with the fast drop in density at small radii for collapsed and virialised structures which results in 
a strongly suppressed contribution to the volume integral of the gradient squared. 

Therefore, we find it to be likely that in coupled quintessence scenarios the clustered matter 
distribution on supercluster scales has a relatively large impact on the average equation of state 
of dark energy, provided the interaction is of at least gravitational strength. 

It is also important to note that the results of our investigation can be taken to be conser- 
vative for the following reasons. As a simplifying assumption, we have neglected the anisotropic 
nature of structure formation by assuming a spherically symmetric matter density profile for the 
superclusters. However, the shape of structures still in the process of formation is triaxial and 
only becomes closer to spherical in the process of collapse and relaxation. Adding to this, while 
structure is infalling, significant internal substructure is still recognisable. Accordingly, within the 
virialised constituents of the supercluster, the scalar field value can be expected to decrease, while 
outside it relaxes again. Within the binding radius of the supercluster, therefore, the coupled 
scalar field will exhibit an oscillatory behaviour. Since the average scalar field equation of state 
is, importantly, only sensitive to the square of the gradient, both positive and negative scalar field 
gradients will act in the same way, to increase the average equation of state. Further study in this 
direction is underway [63]. 

5.4 Comment on the Strength of the Scalar/Matter Coupling 

In accordance with previous numerical studies in the literature [56, 57, 58], we have demonstrated 
analytically that in the present universe the scalar field evolution is non-adiabatic within com- 
pact structures ranging from galactic up to supercluster scales. We therefore found, for coupled 
quintessence models in general, that the relative importance of the average gradient squared of the 
scalar field compared to the average potential grows quadratically with the coupling strength p. 
While this result is largely independent of the scalar field potential V{(j)), possible effects on the 
average equation of state of the quintessence field thus strongly depend on the coupling strength. 
In this section we comment on the implications for various coupled dark energy scenarios of the 
model-dependent constraints on /3 that result from current cosmological observations. 

These bounds strongly depend on the dynamical evolution of the field, controlled by the 
interplay between the scalar potential - and, in particular, its steepness - and the form and 
strength of the scalar/matter coupling governing the effective potential V^s. In the absence of 
a global minimum of T4ff, the ultra-light mass of the scalar field is independent of the local 
matter density, leading to a long-range fifth force between the coupled matter species on all scales. 

In this case, a coupling between cold dark matter and the scalar field is restricted to 5% — 20% 
of gravitational strength by current observations [20, 21, 22]. Therefore, our results demonstrate 
that, even on supercluster scales, the average gradient squared of the scalar field is at least two 
orders of magnitude lower than the average scalar field potential. Within these coupled scenarios, 
deviations from the homogeneous equation of state arc then negligible. 

However, in the case when l^g allows for a minimum, by virtue of the chameleon mecha- 
nism it is possible that in high density regions the field acquires a large mass such that the 
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range of the fifth force becomes undetectably short [25, 26, 27]. Stringent constraints from 
local tests of the equivalence principle as well as the solar system can thereby be evaded even for 
couplings /3 ^ 1 [28, 29, 31]. Even so, on cosmological scales it is possible that the scalar field 
has a low enough mass to affect the non-linear formation of matter structures in a complex man- 
ner [56, 57, 58]. More precisely, the scalar field acts as a potential for the fifth force, which itself 
depends on the underlying matter density in a different (non-linear) manner to the gravitational 
potential^. Taking into account the spatial distribution of the scalar field for such a scenario, 
numerical N-body simulations have been performed to determine the impact of a scalar/matter 
coupling (3 of gravitational strength on the power spectrum [56, 57, 58], the mass function [57, 58] 
and the internal profiles of dark matter halos [58]. One of the main conclusions is that on linear 
and highly non-linear scales, these strongly-coupled chameleon- like models are virtually indistin- 
guishable from ACDM [57]. However, as a result of the enhanced growth in structure due to the 
fifth force, the matter power spectrum was found to be significantly increased with respect to 
ACDM on intermediate scales relevant for galaxy clusters. The authors of Ref. [57] conclude that 
future surveys operating on scales k = 0.1-10 /iMpc~^ will be needed if we are to further tighten 
constraints on the coupling below gravitational strength. 

To summarise, in rather weakly coupled dark energy scenarios with /3 <C 1, the effect of the 
scalar matter coupling on the average equation of state of the scalar field is negligible, as could be 
expected. However, this is not the case in chamclcon-like cosmologies which allow for couplings of 
gravitational strength between the scalar field and dark matter or all kinds of matter. As we have 
demonstrated in the last section, deviations on the level of a few percent from the homogeneous 
equation of state are conceivable. 

In the next section, we will provide an outlook for the implications of our results for even less 
constrained theories such as coupled neutrino quintessence [32, 33, 34, 35, 36] with scalar neutrino 
couplings ^ 2> 1. 

6. Outlook — Couplings of Super-Gravitational Strength 

In this section we briefly consider the implications of our results for dynamical dark energy scenarios 
in which a quintessence field couples to some dark species with a coupling strength much stronger 
than gravity, /3 S> 1, such as in coupled neutrino quintessence [32, 33, 34, 35, 36]. 

In stark contrast to a ACDM cosmology, the formation of structure in the coupled species 
does not cease in the era of dark energy domination. This is due to the fact that gravity, and 
therefore the time scale for gravitational collapse, is directly connected with the evolution of the 
scale factor. However, the attractive force mediated by the scalar field between the coupled matter 
is independent of the background cosmology. Consequently, it allows for structures forming up to 
the present time. As a direct consequence, the characteristic size of these lumps can in principle 
be much larger than superclustcr scales, 80/i^^Mpc ^ Tobj ^ O.IH^^ (depending on the specific 
model assumptions) [50]. Accordingly, even only a few of these objects can take up a very large 
fraction of the total Hubble volume Vh- 



simple rescaling of the gravitational constant G has been shown in certain regimes not to be physical [56, 57, 58]. 
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Figure 5: w as a function of /3 for a volume of 3% (left) and 5% (right) taken by bound structures 
with radius rvir = 20h~^ Mpc and a mass of A/ = 5 x IO^^Mq (solid), M = 10^^ Mq (long-dashed) and 
M = 2 X lO^^^Afo (short-dashed). 

Clearly, as in our current study, the specific functional form of the spherically averaged density 
profile will have a large impact on the size of possible corrections to the average equation of state 
of the coupled scalar field. The authors of [50] found approximately a high-concentration NFW 
profile for neutrino lumps forming in a model with a high coupling strength, |/?| =52 and current 
neutrino masses of 2 eV. 

For highly concentrated structures, the normalised gradient squared Q resulting from Eq. (4.15) 
without employing the relations between c and in Eq. (4.4) is 



To get a crude estimate for the impact of such structures on the cosmic equation of state, we plot 
in Fig. 5 the resulting w as a function of /? S> 1 for typical class of high mass objects with an NFW 
like matter profile. 

Even assuming an occupancy of only 3%, we can see that a strength /3 ~ 50 will produce an 
average equation of state —0.99 < w < —0.92, depending on the ratio Mvir/r^jj,; for an occupancy 
of 5% the average equation of state would be —0.99 —0.88. Therefore, compared to the case 

of couplings of (sub-)gravitational strength, much larger deviations to the homogeneous equation 
of state seem feasible in these kinds of scenarios. 

In the case of an isothermal matter distribution the result can be expected to be yet larger, 
but it will again strongly depend on the model-dependent mass and radius of the object, as well 
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as on its concentration. In this case it should also be verified whether the increase in /3 allows for 
an adiabatic field evolution. 

The possibility of finding notable deviations from the homogeneous case clearly suggest that 
this issue warrants further study. In particular, this might imply a source of new observational 
constraints on the coupling strength. 

7. Conclusions 

In this paper, we have investigated the impact of a coupling /3 between a quintessence field and 
clustered matter on the average scalar field equation of state across the universe w. For this 
purpose, we have derived its general form which turns out to crucially depend on the ratio of 
the average gradient squared to the average potential of the field. In the static limit, this ratio 
vanishes for a perfectly homogeneous uncoupled scalar field, implying To ~ — 1. However, this is 
not the case in the presence of sizeable spatial variations in the field, where in the static limit the 
equation of state could reach an as3miptotic value of uJ 2± —1/3. 

As a first step, we considered isolated bound objects with non-trivial density profiles character- 
istic of their size and evolutionary state. Solving the equations of motion for the field analytically 
for a variety of field potentials, we determined the spatial profile of a quintessence field induced by 
the coupling to matter. To model objects on galactic and cluster scales, we assumed a matter pro- 
file of the NFW form, while on supercluster scales we considered an isothermal matter distribution 
oc modified at small and large radii. 

A first crucial result is that, within bound objects on all scales, the absolute change in the 
value of the coupled scalar field is much smaller than the fundamental field value which is of order 
Mpi. An immediate consequence is that, for a typical quintessence potential, the average of the 
scalar field potential over the volume of any bound object is largely independent of the strength of 
the coupling. Furthermore, this implies that the correction to the uncoupled value of the potential 
is negligible. 

However, the average gradient squared of the field grows quadratically with the strength of the 
coupling. Importantly, its magnitude furthermore exhibits a strong dependence on the functional 
form of the matter distribution; in the case of the NFW profile, the corresponding volume integral 
is strongly suppressed both at small and at large radii. As a result, the average gradient squared 
is at least a factor 0(10^)/?^ smaller than the scalar field potential, even for the largest and most 
massive galaxy clusters. 

More significant results can be found in the case of galaxy superclusters. Fully evolved super- 
clusters obey an isothermal density profile across much of their radii, proportional to r~^. The 
average gradient squared consequently receives almost equal contributions at all radii and, as a 
result, for typical parameters it can be at least j3'^ times larger than the scalar field potential. 

This result depends crucially on the shape of the profile; however, unlike fully virialised objects 
such as galaxies and galaxy clusters, superclusters are still forming. Their spherically averaged 
density profile then depends on their evolutionary state and will deviate from the isothermal 
at small and large radii. To probe the dependence of the results on the evolutionary state of 
the superclusters we considered two representative test objects from simulations, one relatively 
unevolved and one much more so. 
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In both cases, we found the matter in the core region to be well-fitted by a Gaussian distri- 
bution, while at large radii the slope of the distribution dropped faster than r~^. We built an 
analytical model for the density profile, stitching a Gaussian core to an extended isothermal tail. 
We could then solve the equations of motion both analytically and numerically, allowing the fol- 
lowing general conclusions. At small and medium radii, the contributions to the average gradient 
squared of the coupled field resulting from the analytical fit agree very well with the findings from 
the numerical data. However, as soon as the matter profile drops faster than r~^, the integrand of 
the volume integral is no longer approximately constant, but instead decays. Accordingly, further 
evolved superclusters result in larger average field gradients, since the isothermal profile provides 
a better fit for a larger range in radius; in our particular cases the average gradients squared 
were approximately an order of magnitiidc greater for the more evolved object than for the less 
evolved. Adding to this, for a power-law density profile, the magnitude of the average field gradient 
squared grows quadratically with the supercluster mass, but is suppressed by the fourth power of 
the radius. 

Taken across the Hubble volume, these results suggest that the greatest corrections to the 
homogeneous equation of state in coupled quintessence scenarios are expected to originate from 
further evolved, medium sized superclusters; these objects are both large and numerous enough 
to fill a reasonable fraction of the Hubble volume while additionally inducing relatively large de- 
viations from the homogeneous behaviour. We could therefore infer from our supercluster case 
studies that a scalar / matter coupling of gravitational strength allows for an increase in the average 
equation of state of the universe of up to ~ 3%. Let us emphasise that this is a conservative esti- 
mate, since substructures will presumably induce oscillations in the field gradients; both positive 
and negative gradients will increase the result. 

In coupled dark energy scenarios which allow for couplings of gravitational strength, the 
global equation of state of the scalar field across the Hubble volume will therefore be displaced 
from W ^ —1. As it is this equation of state that one would expect to employ in the Priedmann 
equations, the global acceleration rate in coupled quintessence models will therefore be less than 
would be naively expected. 

It should also be noted that particular coupled quintessence scenarios which allow for super- 
gravitational strength couplings - to neutrinos, for example - can produce even larger deviations, 
since the average gradient squared grows quadratically with the coupling strength. We briefly 
considered this scenario using the NFW profile and found even for this distribution that corrections 
across the entire Hubble volume of up to ~ 10% seem feasible for typical model parameters. 
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1. Introduction 

The discovery of cosmic acceleration requires a source of new physics which may come as a modifica- 
tion to gravity, but is entirely consistent with a cosmological constant or quintessence, a dynamical 
dark energy [1, 2, 3].^ Current observations suggest that dark energy comprises approximately 
70% of the energy content of the universe [7, 8] . 

The second major ingredient of modern cosmology is the dark matter, making up approxi- 
mately 25% of the universe's mass [7, 8]. Yet, remarkably, the existence of both of these dominant 

^ Other suggested alternatives include the backreaction of inhomogeneities on the global expansion rate [4]; recent 
implementations of this have individually produced a small dark energy [5] and a large variance in the Hubble rate 
[6] but not provided a unique dark energy source. 
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components of the universe is only deduced from their gravitational interactions. While a funda- 
mental microscopic picture is still lacking, a unified description of the dark sector might exist. 

One way to test such a proposition is to explore the phenomenological consequences of possible 
further non-gravitational interactions between a dynamical dark energy component, an ultra-light 
quintessence field, and dark matter [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]. As an immediate 
consequence of the interaction, the scalar field mediates a long-range force between the dark matter 
particles. As a result, for some scenarios strong constraints on the strength f3 of this dark force have 
been derived from dark matter dynamics [20, 21] and from the modified evolution of cosmological 
perturbations and subsequent impact on the CMB [22] , constraining the coupling to 5% — 20% of 
gravitational strength. 

On the other hand, there arc several observational motivations for dark interactions beyond 
gravity which challenge the ACDM concordance model. These include anomalously large bulk 
flows on scales of 100/i~^Mpc and beyond [23], and higher supercluster densities [24] and sparser 
voids [18] than expected. 

Certain types of adiabatic coupled dark energy scenarios [25, 26, 27] have been shown to 
naturally allow for strong couplings /3 S> 1 [28, 29, 30, 31] between a light scalar and all forms of 
matter - even baryons - by means of the so-called thin-shell or chameleon mechanism. This applies 
if the field is kept in a minimum of its effective potential inside an overdense region and implies 
that deviations from the gravitational force law arc shielded in the interior of bound objects. 

Another alternative for scenarios with couplings beyond gravity recently studied in the liter- 
ature is coupled neutrino quintessence [32, 33, 34, 35, 36]. The scalar /neutrino interaction results 
in growing neutrino masses which can trigger the onset of cosmic acceleration [36] . 

Since all viable dark energy models must yield extremely similar background behaviour, there 
has been intense study into alternative methods to discriminate between individual homogeneous 
and inhomogeneous models at a perturbative level. An uncoupled quintessence field with negligible 
self-interactions is not sensitive to the presence of inhomogeneities in the universe, and thus remains 
almost perfectly homogeneous. However, within overdensities the presence of a scalar/matter 
coupling causes the field to collapse along with the dark matter which itself feels an additional fifth 
force. The resulting modified growth of perturbations has been analysed in the linear [37, 38, 39, 40] 
and in the non-linear [41, 42, 43, 44, 45, 46, 47, 48, 49, 50] regimes. 

However, during structure formation, the dark energy provides a subdominant component to 
the energy budget of the universe, thereby limiting the size of the predicted effects. At late times, 
not only does the dark energy dominate the universe, but the matter distribution is significantly 
non-linear. It then seems natural to expect further sizeable effects at a non-perturbative level, 
emerging from non-standard scalar/matter interactions in the recent universe. In this paper, we 
therefore focus on possible corrections to the homogeneous equation of state of a coupled scalar 
field resulting from its non-trivial spatial profile acquired within clumped matter. 

As a result of its formation and collapse, a gravitationally bound structure in the universe 
is characterised by a non-trivial density profile with a shape which reflects its evolutionary state. 
By the current epoch, galaxies and galaxy clusters have had time to develop into fully virialised 
structures and, according to both observations and simulations, are well-described by a Navarro- 
Frcnk-White density profile. Even though these objects have a sizeable abundance, their total 
volume makes up a relatively small fraction of the Hubble volume. While it is certainly conceivable 
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that local fluctuations in the coupled quintessence field could produce noticeable corrections to 
the homogeneous equation of state from such objects, we might expect the greatest corrections 
to result from the largest and most massive known structures in the universe, the superclusters. 
However, in contrast to fully relaxed objects on smaller scales, superclusters are still currently in 
the process of forming. This implies that the spherically averaged density profile of a supercluster 
depends on its evolutionary state. In the future, it will tend for all superclusters towards an 
isothermal form proportional to at all radii [51]. However, at present the matter profiles of 
less evolved superclusters are typically shallower at small radii, but exhibit a steeper slope at large 
radii [52, 51]. 

Our first aim in this paper is to analytically solve the scalar field equations for different kinds of 
field potential and sourced by matter comprising of isolated, compact objects with density profiles 
of Navarro-Frcnk- White and of isothermal form, corresponding to galaxies and clusters, and to 
superclusters respectively. The resulting spatial profiles of the field will allow us to determine the 
average equation of state of the scalar field within such bound structures. 

Our second aim is to carefully consider the impact of the evolutionary state of superclusters 
on the average equation of state of a coupled quintessence field. For this purpose, we will employ 
the results of numerical simulations to consider two case studies in different stages of formation. 
We will compare our numerical findings with those for an analytical fit to the simulated data, 
namely a Gaussian profile in the core stitched to an extended isothermal tail. 

Taken across the Hubble volume, these results will allow us to estimate the overall effect on 
the global equation of state resulting from the induced inhomogeneous nature of the coupled scalar 
field. 

After setting the stage for our analysis in §2, we define the average equation of state of an 
inhomogeneous scalar field. In §3.1 we provide a model-independent qualitative discussion of 
the scalar field dynamics in the presence of a coupling to clustered matter. In §4 and §5.1 we 
analytically determine the spatial profile of a scalar field coupled to bound matter structures 
with an NFW and an isothermal one profile, respectively, and calculate its average equation 
of state within such objects. In §5.2 we construct an analytical model for simulated density 
profiles of representative superclusters in different evolutionary states. In §5.3 we will discuss and 
compare our analytical and numerical results and estimate the implications on the equation of 
state averaged across the Hubble volume. §6 contains an application of our results to coupled 
quintessence scenarios with interactions beyond gravity and in §7 we conclude. 

2. Setting the Stage 

Our aim is to investigate the dependance of the average equation of state of a quintessence field ^ 
on its coupling to clustered matter. 

2.1 The Matter Distribution 

In the literature, various viable dark energy scenarios have been considered in which a light scalar 
field couples either to dark matter [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19], to neutrinos [32, 33, 34, 
35, 36] or democratically to all forms of matter [25, 26, 27]. In general, all of these types of coupled 
matter have by the present time developed into bound compact objects and were shown to exhibit 
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similar characteristic density profiles [46, 53, 50]. In order to deduce results valid for different kinds 
of coupled quintessence scenarios, we will therefore model the matter distribution as individual, 
isolated, spherically-symmetric clumps of a variety of sizes embedded in a background vacuum; we 
will consider bound objects of galactic size, as well as clusters of galaxies and superclusters. This 
is in contrast to the usual approach where overdensities are embedded in a constant background 
density. This implies that we will neglect non-linear interactions between the objects; since the 
typical distance between each of these bound objects is typically significantly larger than their 
size, this can be justified first approximation. 

We are likewise taking spherical symmetry to hold for objects across a wide range of scales; for 
virialised objects up to cluster scales we will assume a Navarro- Prenk- White (NFW) profile [54], 
while for superclusters we investigate the characteristic isothermal profile along with modifications 
to it at small and large radii (see [52, 51] and references therein). Spherical symmetry is a good 
approximation on smaller scales where objects are fully virialised. When averaging across a large 
number of objects, we can assume it can also be justified on larger scales, and we would expect 
our calculations to provide reliable order-of-magnitude estimates. 

For the purposes of this paper it suffices to approximate the geometry of space as Minkowski, 
Qnv = iJni, = diag(— 1, 1, 1, 1), implying that both the Newtonian potential and the backreaction 
due to the energy density in cf) stay small everywhere. Even though wc will be taking averages 
across the Hubble volume, this approximation is valid since we take the bound objects to be 
isolated from one-another; each resulting average is therefore taken only across a patch of space 
small compared to the Hubble volume. 

In addition, we assume the field dynamics to be static, since the characteristic time scale for 
a dark energy candidate is of the order of the Hubble time. Furthermore, in order to probe the 
dependence of the local profile of cf) on the fundamental potential of the scalar field, we will consider 
two fiducial examples, an inverse-power law and an exponential potential. 

2.2 The Average Equation of State 

In this section, we isolate the relevant quantities that define the average equation of state of a 
scalar field with a self- interaction potential V {(/)). Since an interaction with clustered matter 
will induce local variations in the field value, we do not a priori assume that spatial gradients 
vanish under averaging as they do in the usual homogeneous approximation. 

When considering an "average equation of state", there are three definitions we can readily 
imagine. The equation of state of a scalar field is defined from its pressure and energy density, 
which themselves depend directly on the scalar field itself. We can therefore construct an average 
equation of state wi by forming the local equation of state and averaging this across a spatial 
volume; or we can average the scalar field ^ itself and define the average equation of state W2 from 
the general forms of the pressure and energy density; or we can average the pressure and energy 
density independently and define an average equation of state W3 from these. The three definitions 
can be written 




(2.1) 
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and they coincide for a homogeneous field. Since the relevant thermodynamic quantities are the 
pressure and the energy density rather than the scalar field or the local equation of state, we choose 
the last approach. 

Varying the scalar field Lagrangian with respect to the metric and averaging over the Hubble 
volume Vh, the average pressure and energy density of the field can be defined as 



P4> 



d^x 



2^ 



-A2 . 



+ 2 i^^r + ^i^) 



l-TT 1 



In the static limit ^ 2± 0, the average equation of state w therefore takes the form 



w 



where 



_ Pi 
Qct> 



-i-H_g 



(2.2) 



(2.3) 



(2.4) 



Accordingly, we can already isolate two exact limiting cases: for Q S> 1 the scalar field dynamics 
will be gradient dominated, leading to w — )• —1/3; conversely, for Q ^ 1 the field gradient is 
subdominant compared to the potential, and w ^ —1. This implies, therefore, that in the static 
limit -1 < ItJ < -1/3. 

Gradients of the scalar field originating from an interaction with clustered matter can thus, in 
principle, lead to sizeable corrections to w as compared to its uncoupled, homogeneous analogue 

Wh, 



_ P0 
^h = — 7=r 

Q<l> W 



-V{<f>) 



(2.5) 



where again the last equality applies in the static limit. 

To transfer the calculations for the equation of state of the coupled field within a bound object 
to an average across the Hubble volume it is necessary to make some assumptions concerning the 
distribution of matter sources. We will separate the bound objects into distinct classes, each of 
which produces a particular scalar field profile and therefore the same equation of state. If i 
denotes a particular class of objects with radius r-Qbj,?, then Eq. (2.3) can be approximated by 



W 2±= 



(2.6) 



with 



7^,• = A^, 



Vobj,i 



Qi 



^ob j , i 

_o 

''"obj ,i 

/ drrV((/)(r)) 





(2.7) 
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Here Ni is the number of objects in the class i, we have assumed Q — in the intergalactic medium, 
and angle brackets refer to quantities averaged over the volume of an isolated object Vobj.i- Note 
that TZi is a measure for the volume dilution of the contribution Qi to the sum. We can then 
immediately deduce that notable differences between Eq. (2.7) and Eq. (2.5) could well emerge 
from a class of objects for which the gradient contribution Qi is large enough to compensate for 
the volume dilution implied in TZi. 

3. Dynamics of the Scalar Field 
3.1 The Equation of Motion 

As is clear from the previous section, to estimate w we must determine the spatial profile of the 
scalar field (f> resulting from its interaction with clustered matter. While we assume standard dark 

matter density distributions for the coupled matter species, wc will consider different possible 
forms of coupling functions leading to distinct dynamics of the field. By this means our results can 
be generalised to hold for the various different kinds of coupled scalar field scenarios mentioned in 
§2. 

In general, in the presence of an interaction of strength /3a between ^ and a > 1 non-relativistic 
matter species, the dynamics of (f) are governed by an effective potential of the form 

yes{^,r) = V{<j>,r) + ^Pa{cP,r), (3.1) 

a 

where the first term arises from self-interactions of the quintessence field, generating a monotonically- 
decreasing potential ^(0) of a runaway form. The second contribution to V^ff originates from the 
<p dependance of the masses of the coupled species, ma{(f>), introduced by the interaction. Since 
the corresponding energy density pa is a function of the particle mass, as a direct consequence 
Pai'iTT'ai^')) Contributes to the effective potential for 0. 

For simplicity, we assume a universal coupling of strength P = Pa and an effective energy 
density p leading to an effective potential of the form 

yeff(0,r) = F(</.,r)+p(0,r), (3.2) 

which encodes the effects of one or more coupled species on the scalar field potential. 

Note that depending on the sign of the coupling function dlog(m((;l!)))/d0, we must in general 
discriminate between two cases which lead to distinct dynamics. Namely, a positive sign allows 
the effective potential to exhibit a minimum, while a negative sign implies a (monotonic) decrease 
of VeS with (p. 

We take for simplicity the mass of the coupled matter to depend only very weakly on changes 
in the scalar field value such that it can be taken to be approximately constant even within 
environments with strongly varying matter density - that is, p{(t>,r) 2± p{r). Under the further 
assumption of spherical symmetry for the matter profile, the dynamics of the field in the static 
limit are described by the equation of motion 

0"(r) + -0'(r) = V^{4>, r) + /3p(r) = V,s,^{ct>, r) , (3.3) 
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where here and in the foUowing primes denote derivatives with respect to r, a subscript , </> is a 
derivative with respect to (j), and the last equahty results from Eq. (3.2) and a coupling function 
of exponential form. 

It will transpire that our results are not sensitive to the form of the potential V{(f)). We 
demonstrate this by employing two fiducial examples. The first is of exponential form, 

V{(t>) = M^^e"^^^ , (3.4) 

where Mpi denotes the reduced Planck mass and a > 10 to fit observation [55]. 
Our second fiducial potential is an inverse power law potential of the form 

n<t>) = , (3.5) 

where n is a positive constant and M an intrinsic mass scale. It fixes the behaviour of in the 
intergalactic medium and has to be compared to the scale of dark energy in a homogenous universe. 
However, unless specified the potential should be taken to be general. 

3.2 Qualitative Description of the Solutions 

Let us start our analysis with a model-independent description of different possible kinds of ap- 
proximate analytical solutions to Eq. (3.3). 

We first describe the initial conditions that we assume to hold for any effective potential. 
Since the field equation in Eq. (3.3) is a second order differential equation, we have to specify two 
boundary conditions for 0. Firstly, to guarantee a non-singular behaviour of ^ at the origin, we 
require as usual 

^' = at r = . (3.6) 

To set the second boundary condition, we note that in our setting the bound matter clumps are 
embedded in a vacuum, and so the matter density outside of these asymptotes to zero and not to a 
non-zero background value as it would for an idealised homogeneous universe. We must accordingly 
take the field to be essentially free in the intergalactic medium, that is = ^free- Therefore, in 
'empty space' its dynamics are solely governed by its self-interaction potential T^(0free)- Our second 
boundary condition is then 

^ — >■ 0free as r — >■ oo . (3-7) 

To a good approximation, the solutions for ^ turn out to be independent of the scalar field 
potential ^ ((/»). However, there is one situation in which this is not the case, where (j) — ^min? 
with (/>min being the field value at the minimum. This can occur when the following conditions are 
satisfied: 

VeS,4> = 0, such that V^s exhibits a minimum (3.8) 
0i — 0min(O) <C ^min(O), where (pi is the value of (f) at the origin (3.9) 

^ mri){r), which we refer to as the adiabaticity condition. (3.10) 



If these conditions hold, then the equations of motions allow for an adiabatic solution to Eq. (3.3), 
valid in the range < r < rtick for rkick < '"obj; where the field instantaneously adjusts to spatial 
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changes of the minimum Vcg^^. Here, rtick denotes the radial coordinate at which the fields gets 
kicked out of the minimum. The adiabaticity condition corresponds to the requirement that the 
matter density changes on scales much larger than the Compton wavelength m^^ of the scalar 
field in that region. 

Since in this case the source term is initially T4ff,<^ 0, and the mass of the field inside the 
object is large, m^(0) S> mfree; perturbations in the field will be exponentially suppressed. As a 
result, if r < rtick and Eq. (3.10) is valid, the matter source term p(r) can be treated as a constant 
in the scalar field equation of motion [26]. 

As soon as r > rkickj condition (3.10) is violated and the field gets displaced from the minimum, 
and the contribution from V^f, to the field equation, Eq. (3.3), becomes negligible with respect to 
the matter density p, 

^->V.. (3.11) 

Note that if the condition (3.9) is violated, this implies that the field is already sufficiently 
displaced from the minimum at the origin such that rkick = and its evolution is non-adiabatic 
from the start. 

As it will turn out, for the density profiles and potentials discussed in this paper, Eq. (3.10) is 
already violated for very small r. Accordingly, the field is immediately kicked out of the minimum 
of VeS even in scenarios that in principle might obey all three conditions. As a result, the evolution 
is effectively non-adiabatic across virtually the entire object. 

In general, therefore, in the scenarios we will consider we are lead to the second type of 
approximate solution to Eq. (3.3). This solution is valid as long as Eq. (3.11) is fulfilled and 
implies that Eq. (3.3) reduces to 

r ^!^. (3.12) 

r Mpi 

Following the approach described in [28], we have verified that this approximation stays valid for 
< r < Tobj in the scenarios under investigation in this paper. 

Note for later reference that Eqs. (3.6) and (3.12) yield the following result for the gradient 
squared of (f), 








independent of the form of the scalar potential V((p). 

Averaged over the volume Vobj = 47rrQj^j/3 of the object, this gives 



2 




((V<^)^) = AttV-' ( ) / drr-' I / drTV(r) I . (3.14) 



Finally, outside of the object the dynamics of the field are only governed by its self-interaction 
potential V with <C 1. Accordingly, in this regime the field equation reduces to 

</." + -0'~O, (3.15) 
r 



- 8 - 



which gives the solution 



(7g(-"*free(»'-»-obj)) 

(t){r) = (t)hee • (3.16) 



Here we have employed that the Compton wavelength of the field in the intergalactic medium, 
m^gg, is much larger than the radius, Tobj, of the object. Furthermore, the constant C can be 
determined by requiring continuity of (f) and d^/dr at r = robj- 



4. NFW matter density profile 

We now turn to objects of galactic and cluster scales, for which we employ an NFW matter density 
profile [54] , which is a phenomenological fit to the density distribution of dark matter halos in CDM 
cosmological simulations, 

with 

c^pc ^ rvir . 

3 log(l + c) - c/(l + c) ' r, ' ^ ' 

and we have defined the dimensionless variable R = r/rs, where < R < c, with r being the radial 
coordinate and the characteristic scale of the halo. Together with its virial radius rvir, ?'s defines 
the concentration parameter c. Furthermore, the characteristic density ps is proportional to the 
the critical density of the universe, pc, and to the virial over density, A, which can be computed 
as described in [56]. While the two parameters ps and are in principle independent, they where 
shown in simulations [57] to be correlated through the concentration parameter c, the critical 
density pc and the virial mass Myir, 



Mvir = 4Trpsr, 



log(l + c) - 



1 + c 



(4.3) 



Note that typical virial masses of galaxies and galaxy clusters are Myir — lO^'^ — IO^^Mq and 
Mvir ~ 10^*^ - W^^Mq, respectively. 

In the following, we will employ the relations 



13.6 f J:^^ V°''' , r.^8.8(-^] kpc. (4.4) 



As discussed in the last section, the behaviour of the coupled scalar field inside a matter lump 
crucially depends on the magnitude of its interaction range m^^ compared to the rate of change 
of the matter density, see Eq. (3.10). 

In the following subsection we derive the range of validity of the adiabatic approximation to 
the field equation for the two fiducial potentials in Eq. (3.4) and Eq. (3.5). In addition, we assume 
the matter coupling function to be an increasing function of the field value such that the effective 
potential exhibits a minimum. 
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4.1 Adiabaticity Condition for an NFW Matter Profile 

For the NFW distribution the adiabaticity condition (3.10) becomes 



p'{R) 



p{R) 



^ —R ^ , (4.5) 



Vs R{1 + R) rs 



where the last equality assumes <C 1. As a result of the cusped structure of the matter density, 
Eq. (4.5) diverges at the origin. This already suggests that the field will be immediately displaced 
from the minimum, as we now explicitly verify. 

In the case of the exponential potential, we find small excitations about the minimum of the 
potential to be given by 

(a/3p. 



m*W = «(l + ii)-'iJ-4. (4.6) 



The resulting adiabaticity condition therefore takes the form 

1 -1 

1 9 A/fr., \ 2 

(4.7) 



1 

jj 1,1/^ aA(/)NFW 
-ft <c 



3 18 V ^Pi J V aA0NFW / 
with 

^*-W.^.C(10-)Mp,(^)°", (4.8) 

where the last equality results from Eqs. (4.3)-(4.4) and holds for c > 1. As it will turn out later 
on, A0NFW is a measure for the absolute change of the scalar field value within the object, which 
for typical model parameters is much smaller than the characteristic mass scale in the potential, 
A(/)NFw a/Mpi. Eq. (4.7) therefore cannot be solved for any i? > and the field evolution is 
non-adiabatic from the start. 

In the case of the inverse-power law potential, the field mass expanded around the origin takes 
the form 

m,{R) c V^n-i(^)M-MSil) '^^^ R-'^^) . (4.9) 

Again, for any positive n, Eq. (4.9) diverges less fast than Eq. (4.5) at small radii. 
In this case, the adiabatic condition in Eq. (3.10) close to the origin becomes 

i?< (n + l)(^)n~»M^r7"(A0NFw)~*'^'' < 1, (4.10) 

where the last equality holds even for the largest and most massive clusters of galaxies and realistic 

model parameters. 

Since the slope of the distribution in the innermost region of dark matter halos is still under 
debate (see [59] and references therein), we have also verified that, for forms of the NFW profile 
regularised at small radii discussed in the literature [60, 61], the field is kicked out of its minimum 

very close to the origin. 

Wc therefore conclude that Eq. (3.11) can be assumed to hold within objects up to cluster 
scales. The resulting solutions for the field thereby become independent of the functional form of 
the scalar field potential. 
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4.2 Non-Adiabatic Solution for the NFW Matter Profile 
For an NFW profile, we find that the solutions to the field equation in Eq. (3.12) are 

'log(l + R) 1 



4>m. = <j>h:ee — A^nFW 



R 



1 + c 



^out 



^free - A^NFW log(l + c) 



(4.11) 
(4.12) 



l + cj R 

tree'vir ^ 1- Furthermore, as described in §3.2, we have deter- 



where the result again assumes mfree?^v 
mined the integration constants by matching the inner solution to the outer solution ^out) 
R = c and required continuity of the first derivative of the field at i? = c. 

As we mentioned in the last section, the absolute change in the scalar field value, (f){c) — (p{0) ~ 
A(/)NFW defined in Eq. (4.8) within a compact object with virial radius Tvir, is a measure for the 
total mass Mvir defined in Eq. (4.3), which is interacting with with coupling strength /3. 

Using Eq. (4.11), the volume average of the exponential potential becomes 



{V) ~ Ffr, 



1 + 



NFW 



Mpi 



\ + {l-c) (2-31og2) 



(4.13) 



for c > 1. 

In the case of the inverse power law potential we find 



1 + n 



A^NFW f 1 



+ (l-c) (2-31og2) 



(4.14) 



for c > 1. Note that for both potentials, the matter coupling results in a positive correction term 
with respect to the free field potential T4rce- However, importantly, for typical model parameters 
applying to galaxies and galaxy clusters, it turns out to be negligible, since according to Eq. (4.8) 
A0NFW ^ 0free a-^d A^nfw ^ Mpi/a, respectively. 

For the volume average of the gradient squared we get 



{iVcf>f) = 3c 



-3 f A^NFwV 



dR 



log(l + R) 
R 



1 



1 + R 



3/(c 



/ A^NFW Y 

V crs J 



with 



/(c) = 



1 + c 



log(l + c) 



(4.15) 



(4.16) 



where /(c) ~ O(10~^) for typical values of the concentration parameter c. Note that even though 
the matter density formally diverges at the origin, the gradient squared stays well-defined dX R = 0. 
This is due to fact that the volume integral over the matter density is finite as can be seen from 
Eq. (3.14). 

Accordingly, using Eq. (4.4), the volume average of the gradient squared normalised to the 
free field potential becomes 

\2\ / )ij \ 0.68 



Q 



((V'/')0 

^free 



5.12 X 10-^/3^ 



(4.17) 
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for typical galaxies, and 

for higher mass clusters, such that Q ^ 1 for couplings of (sub-)gravitational strength ,0 < 1 (see, 
however, §6, where wc analyse the results for coupled quintessence scenarios with interactions of 
super-gravitational strength). 

Recall that, according to Eq. (2.7), contributions to W depend on the product of Q and TZ. 
However, for all classes of objects on galactic and cluster scales, the volume suppression amounts 
to at least TZ < O{10^^^)N <^ 1 for typical abundances of these object classes, w therefore receives 
negligible corrections compared to the uncoupled case Wh — —1- 

As a concluding remark, let us note that the small result for Q is a direct consequence of the 
specific functional form of the NFW profile. Namely, according to Eq. (4.15), the integrand of 
the volume average is of 0{R^) for R 1 and of 0{1/R^) for R^ 1. This implies that in both 
regimes the corresponding contributions to the volume average are strongly suppressed and thus 
the average gradient squared only receives sizeable contributions from radii i? ~ 1 <C c. 

However, as we will see in the next section, for the density distribution characteristic for 
objects on supercluster scales, the average field gradient receives almost equal contributions at all 
radii and thus turns out the be much larger. 



5. Superclusters 

In this section we investigate the impact of superclusters, the largest known bound structures in 
the universe, on w in coupled quintessence. Within a ACDM cosmology, superclusters are expected 
to have masses ranging from Mgc ~ lO^^/i~^M0 up to greater than 6 x 10^^/i~^MQ, with sizes 
~ 10/i~^-160/i~^Mpc [52, 62]. Superclusters therefore fill approximately 3% — 5% of the current 
Hubble volume [52]. 

Since structure forms hierarchically, the superclusters we observe in the present universe are 
still in the process of formation, while galaxies and clusters of galaxies have already developed into 
collapsed and virialised structures. Accordingly, depending on their evolutionary state, superclus- 
ters tend to have a much shallower density profile than fully relaxed objects. 

According to simulations, the spherically averaged density profile of superclusters is well-fitted 
by an isothermal profile p{r) oc over a broad range in radius. While the profile is shallower 
than this at low radii, it steepens at larger radii [52, 51]. Since we only have statistical information 
on the exact density distribution, we will proceed in the following manner. 

To establish some intuitive understanding of the behaviour of the field inside a typical super- 
cluster, we will begin by analytically determining the generic profile of a scalar field interacting 
with isothermally distributed matter. This allows us to get a feeling for the magnitude of the 
effects on the global equation of state as a function of the coupling, the mass and the radius of a 
single, standard supercluster. 

Subsequently, we will turn to two simulated supercluster case studies. These allow us to test 
both analytically and numerically the dependence of our results on deviations from the isothermal 
profile for small and large supercluster radii. 
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Finally, we will end our discussion of superclusters by estimating the overall effect on the 
global equation of state. 



5.1 Isothermal Matter Distribution 

Formally, the isothermal profile p{R) oc exhibits a singularity at the origin, where R = r/rsc 
denotes the radial coordinate normalised to the supercluster radius rsc- However, bound objects 
which have not yet relaxed will exhibit a flattening of the matter distribution at some R = A with 
< A <C 1. Accordingly, for A < R < 1 the isothermal distribution expressed in terms of the 
finite matter density p = p{A)A'^ reads 



With the help of Eq. (3.11), the field equation Eq. (3.3) then takes the form 

d'^d) . 2 d</> /3 



di?2 RdR 



P_ 

Mpi i?2 



Assuming both the field value and its first derivative to be continuous at R 
A< R<1 that 



4>m = ^free - A0is 
^free - A0is 



1 - log(i?) - 



A 
R 



R=A/ J 



1-A I 



where 



Mpi 



R=A 



(5.1) 
(5.2) 

A, we get for 

(5.3) 
(5.4) 

(5.5) 



As described in §3.1, we have determined the integration constants by matching the inner solution 
(l)in to the outer solution (j)out, and by requiring continuity of the respective derivatives at i? = 1. 

Note that we can disentangle the different contributions to the average gradient squared by 
investigating the limiting case A0is S> A d(?!>/di?|^^^. Doing so assumes the profile of (j) in the 
isothermal region to be insensitive to its slope in the core region; in general, this slope is associated 
with large uncertainties introduced by the matter profile in the core. Comparing this with the 
analysis in §5.2, where we employ a Gaussian distribution at small radii, will allow us to gain some 
understanding of the influence of the core profile. 

According to Eq. (5.3), for A(^is S> A d(j)/dR\j^^y^ we find the gradient squared of an object 
with isothermal density profile averaged over the volume of the object to be 



{iVct>f) = 3 



IS 



rsc 



1 



/ 



dR(^l 



A 
R 



J AMY 
V ^sc / 



1 + 



21og(^) 
A 



^ 3 (^) ' [1 + 2 log{A)A + 0(^2)1 , 

V ^sc / 



(5.6) 



where the last equality assumes A <^ 1. Note that in this limit the integrand of the volume average 
is of order unity for ^ <C i? < 1 implying that the average gradient squared receives almost equally 



- 13 - 



large contribution at all radii in this range. Contrasting this result with its analogue in the NFW 
case in Eq. (4.15), we observe that the isothermal matter distribution results in an average gradient 
squared which is a factor /^^(c) ~ 0(10^) larger. We therefore arrive at the crucial result that the 
impact of the interaction on the coupled scalar field equation of state strongly depends on the radial 
matter distribution. Since the spherically averaged supercluster profile is typically isothermal over 
a broad range in radius, this result already suggests that the dominant contribution to the average 
gradient squared of the field will originate from this region. 

5.2 Composite Supercluster Profiles 

In the quantitative investigation that follows, we consider two representative superclusters, SC 8 
and SC 98, with masses of Msc ~ 5.4 x lO^^h'^MQ and Msc ~ 3.6 x lO^^/j-^Mo, respectively, 
taken from a simulation performed on the Beowulf Cluster at the University of Groningen [51]. 
While at the current epoch, SC 8 is already a centrally concentrated object, SC 98 is still more 
extended. 

However, both supercluster density profiles exhibit as a characteristic basic shape a high- 
density core with a central density of ~ (10*^ — I0^)pc, embedded within an isothermal power law 
region with slope —2. 

Furthermore, for both of these examples we find that the central region is well-fitted by a 
Gaussian density distribution. Since the two examples represent in some respect two extremes of 
supercluster evolution, it is a reasonable assumption that this is a generic feature. In this section, 
therefore, we model a generic supercluster density profile by stitching together at low radii a 
Gaussian and an isothermal matter density profile. This then allows us to determine analytical 
solutions for the scalar field profile within typical high-mass superclusters which will be compared 
to our corresponding numerical results for the spherically averaged density profiles of SC 8 and 
SC 98. 



5.2.1 Gaussian Matter Distribution 

In the following we determine the field profile (pQ arising from a small Gaussian matter region 
stitched to an extended isothermal tail. The Gaussian and isothermal regions are connected at a 
small radius yl <C 1. 

We parameterise the Gaussian density profile by 

PG(i?)=Pmaxe-^(^^''^% (5.7) 

with 



A 



(5.8) 



We have defined the dimensionless constant Cg = robj/c, a being the half- width of the Gaussian 
distribution with amplitude A. 

In the following subsection we derive the range of validity of the adiabatic approximation to 
the field equation for the two fiducial potentials in Eqs. (3.4-3.5) and will argue that the results 
hold for any realistic analytical scalar field potential. 
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5.2.2 Adiabaticity Condition for a Gaussian Core Profile 

If the effective potential exhibits a minimum and V{(j)) is of exponential form, small perturbations 
about the minimum are described by 



Accordingly, the minimum condition in Eq. (3.10) reads 



v/2 -Iwfii^^) 1 aA(t)G 
R <C — e 2 \2 Mpi J _ ^ I ^ (5.10) 

Cg Cg Y Mpi 

where the last equalities hold for typical parameter values for which aA(^(3/Mpi ^ 1 and Cg S> 1. 
Here W{x) is the Lambert function, which satisfies W{x)e^^^^ = x. 

If the potential exhibits a minimum and V{(l)) is of inverse power law form, we get 

m^{R)=Jc(^]e-^im)(^'^^^\ (5.11) 

with 

K = (nM"+4) 2(n+i) ^ A(t)l}"+'' . (5.12) 
With the help of Eq. (5.7) the minimum condition in Eq. (3.10) becomes 

i^«^e-iMKf$lrK^)^^«l. (5.13) 

Furthermore, the last equalities hold for /C <C 1 and Cg 3> 1 which will turn out to be a good 
approximation for typical parameter values (see §5.3). 

As we will demonstrate explicitly in §5.3, the minimum solution in practice has a very limited 
range of validity, since the evolution turns out the be non-adiabatic very close to the origin. 

We note that these results are in contrast to the findings for compact bodies with top-hat 
profiles which have been employed to analyse solar system bounds arising for chameleon type 
scenarios [25, 26, 28]. The implications of our results, which suggest the absence of a thin-shell for 
objects on galactic up to supercluster scales, will be discussed elsewhere [63]. 

5.2.3 Non-Adiabatic Solution for the Gaussian Core Profile 

In the non-adiabatic case, in the core region, < R < A, the field equation takes the form 

d^ 2#^^ e-^(«c,)^ (5 14) 



with 

where the last equality results from matching Eq. (5.1) and Eq. (5.7) at R = A 



^.e-^^^'^^r, (5.15) 
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Requiring d(f)/dR\j^^Q = 0, with the help of Eq. (5.3) and Eq. (5.5) for < < A we arrive 



at 



(f>G{R) = 0free - ^(f>G 



A<^G = 



max / '"sc 

1 



{ACg 



A(/)is, 



V-l-log(A)A' 



(5.16) 

(5.17) 
(5.18) 



where the relation in Eq. (5.18) guarantees continuity of (p and its derivative at R = A. 

We now have an analytical solution for which allows us to determine the slope of the 
isothermal profile at R = A. 

5.2.4 Isothermal Distributions with a Gaussian Core Region 

In this section, we determine the dependence of {(V(f))'^) in the isothermal region on the slope of 
the field profile in the core region. For this purpose, we require the scalar field profile and its 
derivative to be continuous at the matching point A. Recall that this condition fixes the relation 
in Eq. (5.18). 

Taking the derivative of the field profile in Eq. (5.16) at R = A and inserting it for d^/di?!^^^ 
in Eq. (5.3), we arrive at 



"free 



A(/>is 



1 



log(i?) - |x 



(5.19) 



0out = '?!'free - A^is [1 - Ax] 



R 



where 



1 — erf 



X= 1 + 



ACg_ 

V2 



{Ac, 



-9) 



The scalar field profile in Eq. (5.19) leads to a gradient squared averaged over the volume of the 
object of 

2 



~ 3 



IS 



IS 



i+l^-i)^V + (x^-i)^ 



l+{Uog{A)-^)A + 0{A') 



(5.20) 



where in the last equality we have employed Eq. (5.18) and assumed A <C 1. Noting that for 
X — > 1 we recover the result in Eq. (5.3) with A d0/di?|^^^ <C A^is , we are now in a position to 
estimate the impact of the field derivative at the matching point on the average gradient squared. 
For A <^ 1 and log (A) <C 0, the field gradient at R = A leads to an increase of the result by a 
factor ~ (1 — 2/3\og(A)A) > which amounts to 0(15%) for typical parameter values. 

Accordingly, we find that a small central region with (approximately) Gaussian matter dis- 
tribution has a sizeable effect on the scalar field dynamics. This has to be contrasted with the 
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negligible impact of the corresponding NFW matter distribution at small radii which was found 
to result from its cusped structure. 

For the volume average of the potential in the case of an exponential potential we get 



{V) ~ Ffi 

In the case of an inverse power law we arrive at 

(V) Ffree 



3 / Mpi 



1 + ( ^ - n.^-^^"^ 



3 / ^free 



(5.21) 



(5.22) 



In both cases, the quintessence potential could be expanded around the frcc-ficld potential for 
aA0/Mpi ^ 1 and Acp/cpfiee ^ 1) respectively. As in the case of the NFW profile in §4, we find a 
negligible correction term. 

Note that our result can be generalised to any analytic quintessence potential which allows for 
such an expansion. So long as Acf) is much smaller than the characteristic mass scale of the scalar 
potential, the correction term will remain negligible. 

In summary, for all compact objects from galactic up to supercluster scales, the absolute 
change in the scalar field value induced by the scalar /matter coupling has a negligible effect 
on the average potential of the scalar field. Taken across the Hubble volume, this implies that 
possible corrections to the homogeneous average equation of state can directly be traced back to 
the square of non-trivial field gradients. Note also that as a result of the small present value of 
V 2± Vfree — 10~^^^Mpj, the field gradients only have to be of comparable magnitude in order to 
induce sizeable effects on the average equation of state. 

5.3 Supercluster Case Studies: Analytical and Numerical Results 

In this section we discuss our analytical and numerical results for SC 8 and SC 98, which we take 
as representatives for further and less evolved superclusters respectively. The difference in their 
evolutionary state is reflected in the shape of their spherically averaged density profiles, plotted in 
Fig. 1 together with a Gaussian (for i? <C 1) and an isothermal profile (for R < A) fitted to the 
data. The distribution of SC 8 drops faster than at i? ~ 0.5; in the case of SC 98 the slope 
is already changing at a much smaller radius, R ~ 0.1. Note also that SC 8 has a larger mass 
than SC 98 and therefore its impact on the the scalar field value is larger. This is demonstrated 
in Fig. 2, which shows the scalar field profile for the two superclusters according to Eq. (5.16) and 
Eqs. (5.19-5.20). 

We are now in a position to verify the limited range of validity of the adiabatic solution 
according to Eqs. (5. 10) and (5.13). For our test superclusters we arrive at 

A(/)G = 1.8 X lO-^Mpi (5.23) 

for SC 8, and 

A(j)G = 5.8 X lO^^Mpi (5.24) 

for SC 98. Therefore, for both cases, the adiabaticity condition is already violated at i? ^ 1 such 
that in practise the field evolution is non-adiabatic. 
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In the following, we will compare and discuss our analytical and numerical results for Q as 
defined in Eq. (2.7) and analyse the impact of deviations from the isothermal profile at large radii. 
Finally, based on our results, we will estimate w. 

Taking Fig. 1 as a rough guideline for the relevant scales, according to Eq. (5.19) we find for 
the isothermal matter distribution that 

2 / ~ \ 2 



V, 



free 



30/1-1 Mpcy y42p, 



(5.25) 
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Figure 1: Spherically averaged matter density profiles of the further and less far evolved superclusters 
8 and 98, respectively. Short dashed lines denote the Gaussian core region, long dashed lines the R^^ 
isothermal tail, and solid black lines the numerical results. 



To compare this result with our numerical findings, we insert the fitted radial density distri- 
butions shown in Fig. 1 into the general solution for the average gradient squared in Eq. (3.14)^. 
Fig. 3 shows our numerical and analytical results for the normalised average gradient squared. 



R 



QiR) = ^fr"ee / (^0)' , (5-26) 


and its integrand of as a function of R. As demonstrated in the lower panels of Fig. 3, the final value 
of the normalised average gradient squared for SC 8 is Qs = 1-63, almost an order of magnitude 



■^Note that we have extended the data at small radii with the analytical data from the Gaussian fit in order to 
employ the initial condition in Eq. (3.6) for at R = 
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Figure 2: The scalar field profile (/)(i?)/(/)fioc for the analytie fits for SC 8 (lower curve) and SC 98 (upper 
curve) and 13 — 1. The solutions in the Gaussian region are plotted with short dashed lines, those in the 
isothermal region with long dashed lines, and the external solutions with solid black lines. 



larger than that for SC 98, where Qgs = 0.21. Furthermore, for the former the analytical and 
numerical curves are in very good agreement for all radii, while in the latter case the analytical 
result overestimates Q{R) even at small radii. As can be understood from the plots in the upper 
panel of Fig. 3 of the integrand of Q in Eq. (5.26), the reason is that the isothermal profile overall 
provides a better fit to the profile of the further evolved SC 8. 

According to the analytical result discussed at the end of §5.1, Q receives almost equal con- 
tributions at all radii in the isothermal region, which is reflected by the small slope of the dashed 
curves. As could be expected from Fig. 1, the numerical curves agree with the analytic up to the ra- 
dius at which the numerical profiles drop faster than R^^. However, for larger radii, the numerical 
contributions to the integrand Q{R) show a decaying behaviour. Therefore, to a good approxima- 
tion the corresponding maximal value of Q{R) fixes the total Q, since further contributions are 
negligible. 



Based on our results in this section, we are finally in a position to estimate the overall effect on 
the expansion rate for coupled quintessence scenarios. To achieve this, we note that a particular 
class of high-mass superclusters presumably yields the dominant impact on the expansion rate; 
superclusters tend to become ever less concentrated the larger (and less evolved) they are, as 
reflected by their mean density, 

p = 3- , (5.27) 
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Figure 3: R^{V(t>f (top) and Q{R) (bottom) for SC 8 (left) and SC 98 (right) and /3 = 1. 

which imphes that ((V0)^) oc (/3/Mpi)^ M^^/r^^. Here, the factor of proportionahty depends on the 
radial distribution of mass within the object. Adding to this, the abundance of bound structure 
drops with increasing radius of the object. We therefore conclude that the main source of scalar 
field gradients will stem from medium-sized superclusters which are both large and numerous 
enough that their volume still makes up a reasonable fraction of the Hubble volume (see for 
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comparison the discussion in §2.2 below Eq. (2.7)). 




Figure 4: w as a function of Q for a volume of 3% (left) and 5% (right) and /3 = 1. 



As a result, we assume the total volume taken up by all superclusters (~ (3 — 5)%Vi^) to be 
mainly provided by one representative class of superclusters with a medium size of ~ 30/1-1 Mpc 
and normalised gradient squared Q. According to Eq. (2.7), this corresponds to IZi = TZ ~ 
(3 - 5) X 10^2 

Fig. 4 we plot the resulting w for two different values of 7^ as a function of Q 
together with the corresponding values for w which we get for superclusters similar to SC 8 and 
SC 98. 

From these plots we can see that, if superclusters similar to SC 98 dominate, only small devia- 
tions from IZJ ~ — 1 are to be expected, even should the superclusters fill 5% of the Hubble volume. 
However, for further evolved superclusters of type SC 8, the deviations from the homogeneous 
equation of state are dramatic. For a lower abundance, where the superclusters fill 3% of the 
universe, the equation of state is pushed to ~ —0.98, slowing the cosmic acceleration to some 
degree. For the higher abundance of 5%, the equation of state is driven to ItJ ~ —0.97. Let us 
remark that the authors of [52] have performed a complete analysis of a dark matter simulation 
taking into account the shape, multiplicity, size, and radial structure of superclusters. They found 
that across the entire sample the spherically averaged density profile is roughly isothermal in the 
range 0.3 ^ R ^ 0.7. This investigation suggests that on average, superclusters tend to be much 
further evolved than SC 98 and, in fact, exhibit a similar evolutionary state to SC 8. In light of 
these results, we therefore conclude that an overall correction of > 3% to the homogenous equation 
of state is conceivable in scenarios with scalar /matter couplings of gravitational strength. 
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In summary, as anticipated, in coupled quintessence scenarios the main source for notable 
corrections to w stems from the largest and most massive bound objects known. Since superclus- 
ters are typically less than two orders of magnitude smaller than the Hubble radius, the volume 
suppression is many orders of magnitude lower than that for galaxies and galaxy clusters. 

In addition, the spatial changes in the coupled scalar field value were found to get more 
pronounced for higher masses of bound objects and are also sensitive to its distribution. Namely, 
for superclusters, the dominant isothermal region exhibits a shallower slope than the NFW profile 
discussed in §4 and, importantly, lacks a very highly-concentrated central core (the cusp). Thus, 
according to their evolutionary state, the cumulative mass of an supercluster and its impact on the 
scalar field value change gradually over a much larger range in radius. This has to be contrasted 
with the fast drop in density at small radii for collapsed and virialised structures which results in 
a strongly suppressed contribution to the volume integral of the gradient squared. 

Therefore, we find it to be likely that in coupled quintessence scenarios the clustered matter 
distribution on supercluster scales has a relatively large impact on the average equation of state 
of dark energy, provided the interaction is of at least gravitational strength. 

It is also important to note that the results of our investigation can be taken to be conser- 
vative for the following reasons. As a simplifying assumption, we have neglected the anisotropic 
nature of structure formation by assuming a spherically symmetric matter density profile for the 
superclusters. However, the shape of structures still in the process of formation is triaxial and 
only becomes closer to spherical in the process of collapse and relaxation. Adding to this, while 
structure is infalling, significant internal substructure is still recognisable. Accordingly, within the 
virialised constituents of the supercluster, the scalar field value can be expected to decrease, while 
outside it relaxes again. Within the binding radius of the supercluster, therefore, the coupled 
scalar field will exhibit an oscillatory behaviour. Since the average scalar field equation of state 
is, importantly, only sensitive to the square of the gradient, both positive and negative scalar field 
gradients will act in the same way, to increase the average equation of state. Further study in this 
direction is underway [64]. 

6. Outlook — Couplings of Super-Gravitational Strength 

In this section we briefly consider the implications of our results for dynamical dark energy scenarios 
in which a quintessence field couples to some dark species with a coupling strength much stronger 
than gravity, /3 S> 1. 

In stark contrast to a ACDM cosmology, the formation of structure in the coupled species 
does not cease in the era of dark energy domination. This is due to the fact that gravity, and 
therefore the time scale for gravitational collapse, is directly connected with the evolution of the 
scale factor. However, the attractive force mediated by the scalar field between the coupled matter 
is independent of the background cosmology. Consequently, it allows for structures forming up to 
the present time. As a direct consequence, the characteristic size of these lumps can in principle 
be much larger than supercluster scales, 80/t~^Mpc robj ^ O.IHq ^ (depending on the specific 
model assumptions) [50]. Accordingly, even only a few of these objects can take up a very large 
fraction of the total Hubble volume Vh- 
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Figure 5: w as a function of /3 for a volume of 3% (left) and 5% (right) taken by bound structures 
with radius rvir = 20h~^ Mpc and a mass of Af = 5 x 10^^ Mq (sohd), M = lO^^M© (long-dashed) and 
M = 2 X 10^*^ Mq (short-dashed). 

Clearly, as in our current study, the specific functional form of the spherically averaged density 
profile will have a large impact on the size of possible corrections to the average equation of state 
of the coupled scalar field. The authors of [50] found approximately a high-concentration NFW 
profile for neutrino lumps forming in a model with a high coupling strength, |/3| = 52 and current 
neutrino masses of 2 eV. 

For highly concentrated structures, the normalised gradient squared Q resulting from Eq. (4.15) 
without employing the relations between c and in Eq. (4.4) is 

^ = " '■'%52 J [WH4^) [ lOh-mpc ) ■ ^'-'^ 

To get a crude estimate for the impact of such structures on the cosmic equation of state, we plot 
in Fig. 5 the resulting ItJ as a function of /? ^ 1 for typical class of high mass objects with an NFW 
like matter profile. 

Even assuming an occupancy of only 3%, we can see that a strength /? ~ 50 will produce an 
average equation of state —0.99 —0.92, depending on the ratio Mvir/r^ij.; for an occupancy 

of 5% the average equation of state would be —0.99 ^ < —0.88. Therefore, compared to the case 
of couplings of (sub-)gravitational strength, much larger deviations to the homogeneous equation 
of state seem feasible in these kinds of scenarios. 

In the case of an isothermal matter distribution the result can be expected to be yet larger, 
but it will again strongly depend on the model-dependant mass and radius of the object, as well 
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as on its concentration. In this case it should also be verified whether the increase in (5 allows for 
an adiabatic field evolution. 

The possibility of finding notable deviations from the homogeneous case clearly suggest that 
this issue warrants further study. In particular, this might imply a source of new observational 
constraints on the coupling strength. 

7. Conclusions 

In this paper, we have investigated the impact of a coupling (5 between a quintessence field and 
clustered matter on the average scalar field equation of state across the universe w. For this 
purpose, we have derived its general form which turns out to crucially depend on the ratio of 
the average gradient squared to the average potential of the field. In the static limit, this ratio 
vanishes for a perfectly homogeneous uncoupled scalar field, implying 1(7 ~ —1. However, this is 
not the case in the presence of sizeable spatial variations in the field, where in the static limit the 
equation of state could reach an asymptotic value of uJ 2± —1/3. 

As a first step, we considered isolated bound objects with non-trivial density profiles character- 
istic of their size and evolutionary state. Solving the equations of motion for the field analytically 
for a variety of field potentials, we determined the spatial profile of a quintessence field induced by 
the coupling to matter. To model objects on galactic and cluster scales, we assumed a matter pro- 
file of the NFW form, while on supercluster scales we considered an isothermal matter distribution 
oc modified at small and large radii. 

A first crucial result is that, within bound objects on all scales, the absolute change in the 
value of the coupled scalar field is much smaller than the fundamental field value which is of order 
Mpi. An immediate consequence is that, for a typical quintessence potential, the average of the 
scalar field potential over the volume of any bound object is largely independent of the strength of 
the coupling. Furthermore, this implies that the correction to the uncoupled value of the potential 
is negligible. 

However, the average gradient squared of the field grows quadratically with the strength of the 
coupling. Importantly, its magnitude furthermore exhibits a strong dependence on the functional 
form of the matter distribution; in the case of the NFW profile, the corresponding volume integral 
is strongly suppressed both at small and at large radii. As a result, the average gradient squared 
is at least a factor O(10^)/3^ smaller than the scalar field potential, even for the largest and most 
massive galaxy clusters. 

More significant results can be found in the case of galaxy superclusters. Fully evolved super- 
clusters obey an isothermal density profile across much of their radii, proportional to r~^. The 
average gradient squared consequently receives almost equal contributions at all radii and, as a 
result, for typical parameters it can be at least times larger than the scalar field potential. 

This result depends crucially on the shape of the profile; however, unlike fully virialised objects 
such as galaxies and galaxy clusters, superclusters are still forming. Their spherically averaged 
density profile then depends on their evohitionary state and will deviate from the isothermal 
at small and large radii. To probe the dependence of the results on the evolutionary state of 
the superclusters we considered two representative test objects from simulations, one relatively 
unevolved and one much more so. 
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In both cases, we found the matter in the core region to be well-fitted by a Gaussian distri- 
bution, while at large radii the slope of the distribution dropped faster than r~^. We built an 
analytical model for the density profile, stitching a Gaussian core to an extended isothermal tail. 
We could then solve the equations of motion both analytically and numerically, allowing the fol- 
lowing general conclusions. At small and medium radii, the contributions to the average gradient 
squared of the coupled field resulting from the analytical fit agree very well with the findings from 
the numerical data. However, as soon as the matter profile drops faster than r~^, the integrand of 
the volume integral is no longer approximately constant, but instead decays. Accordingly, further 
evolved super clusters result in larger average field gradients, since the isothermal profile provides 
a better fit for a larger range in radius; in our particular cases the average gradients squared 
were approximately an order of magnitude greater for the more evolved object than for the less 
evolved. Adding to this, for a power-law density profile, the magnitude of the average field gradient 
squared grows quadratically with the supercluster mass, but is suppressed by the fourth power of 
the radius. 

Taken across the Hubble volume, these results suggest that the greatest corrections to the 
homogeneous equation of state in coupled quintessence scenarios are expected to originate from 
further evolved, medium sized super clusters; these objects are both large and numerous enough 
to fill a reasonable fraction of the Hubble volume while additionally inducing relatively large de- 
viations from the homogeneous behaviour. We could therefore infer from our supercluster case 
studies that a scalar /matter coupling of gravitational strength allows for an increase in the average 
equation of state of the universe of up to ~ 3%. Let us emphasise that this is a conservative esti- 
mate, since substructures will presumably induce oscillations in the field gradients; both positive 
and negative gradients will increase the result. 

The global equation of state of the quintessence field across the Hubble volume will therefore 
be displaced from ItJ ~ — 1. As it is this equation of state that one would expect to employ in the 
Friedmann equations, the global acceleration rate in coupled quintessence models will therefore be 
less than would be naively expected. 

It should also be noted that particular coupled quintessence scenarios which allow for super- 
gravitational strength couplings - to neutrinos, for example ~ can produce even larger deviations, 
since the average gradient squared grows quadratically with the coupling strength. We briefly 
considered this scenario using the NFW profile and found even for this distribution that corrections 
across the entire Hubble volume of up to ^ 10% seem feasible for typical model parameters. 
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